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ABSTRACT 


rK 

Variations  in  marine  stratocumulus  cloud  microphysics  during  FIRE  IFO  19S7  are 
observed  and  analyzed  through  the  use  of  NOAA-9/10  AVIIRR  satellite  data  and  air¬ 
craft  in-cloud  measurements.  The  relationships  between  channel  3  reflectance  and  cloud 
microphysical  properties  are  examined  through  model  reflectances  based  on  Mie  theory' 
and  the  delta-Eddington  approximation,  and  reveal  a  channel  3  reflectance  dependence 
on  cloud  droplet  size  distribution.  The  satellite  observations  show  significant  regions 
of  continental  influence  over  the  ocean  through  higher  channel  3  reflectances  resulting 
from  the  injection  of  continental  aerosols  and  the  associated  modification  of  cloud 
droplet  characteristics.  Additionally,  channel  3  reflectance  gradients  across  individual 
cloud  elements  correspond  to  radially  varying  cloud  droplet  size  distributions  within  the 
elements.  Various  mesoscale  and  microscale  features  such  as  ship  stack  diluent  tracks 
and  pollution  sources  are  observed  in  the  data.  Correlations  between  reflectance  values 
and  aircraft  measurements  illustrate  the  potential  of  estimating  cloud  droplet  size  dis¬ 
tribution  and  marine  atmospheric  boundary  layer  aerosol  composition  and  concen¬ 
tration  through  fh'e\ise  of  satellite  data.  Such  an  estimation  technique  may  prove  useful 
in  determining  the  climatic  implications  of  cloud  rellectance  changes  due  to  the  influence 
of  natural  and  man-made  aerosol  sources,  and  provide  a  means  to  assess  the  perform¬ 
ance  of  boundaryTayer  electro-optic  systems.  I  1  r  .  !•  v 
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1.  INTRODUCTION 


A.  ELECTRO-OPTIC  SYSTEM  DEPENDENCE  ON  AEROSOLS  AND  CLOUD 
DROPLETS 

The  molecular  and  aerosol  components  of  the  atmosphere  which  are  within  an  order 
of  magnitude  of  the  operating  wavelength  of  a  military  sensor  or  weapons  system  may 
absorb  emitted  energy,  or  scatter  a  large  percentage  of  it  out  of  the  beam  path,  severely 
reducing  its  efficiency.  Infrared  wavelength  imaging  and  ranging  systems,  currently  de¬ 
ployed  in  numerous  military  applications,  are  shown  by  Milton  et  al.  (1977)  to  be  pro¬ 
foundly  affected  by  variable  distributions  of  atmospheric  aerosols.  Laser  systems  are 
also  strongly  affected,  especially  systems  designed  to  operate  from  ground-based  or 
shipboard  platforms  where  beam  paths  are  confined  to  the  lower  atmosphere.  Accord¬ 
ing  to  Bloembergen  et  al.  (1987),  the  optimum  wavelength  region  for  laser  systems  is 
U.-4-l.O^m.  This  band  avoids  severe  attenuation  from  atmospheric  molecular  compo¬ 
nents  at  short  wavelengths  and  water  vapor  absorption  at  longer  wavelengths.  Atmo¬ 
spheric  aerosols  in  this  band  are  still  strong  contributors  to  energy  losses.  Cordray  et 
al.  (1977)  describes  meteorological  effects  on  a  horizontally  propagating  deuterium 
llouride  chemical  laser  operating  at  3.2/am  in  the  lower  atmosphere,  and  recognizes  that 
aerosol  density,  size  distribution,  and  chemical  and  physical  composition  constitute  most 
of  the  meteorological  parameters  affecting  laser  sensitivity.  The  aerosol  contribution  is 
critical  anu  Cordray  et  al.  recommends  real-time  microscale  measurement  capability  at 
each  sensor  platform  in  order  to  assess  and  predict  the  performance  of  friendly  and 
threat  electro-optic  systems,  and  to  plan  strategy  and  tactics.  A  knowledge  of  the  exact 
aerosol  size  distribution  is  especially  valuable  for  free  electron  lasers  since  they  possess 
a  variable  output  wavelength  which  can  be  adjusted  to  minimize  atmospheric  effects. 
Although  several  high  power  chemical  laser  systems  have  already  been  built  and  their 
use  has  been  under  study  much  longer  than  other  directed  energy  systems.  Bloembergen 
et  al.  (1987)  state  that  these  lasers  are  less  attractive  candidates  for  weapons  systems. 
While  chemical  lasers  are  closest  to  the  required  performance  levels  in  a  relative  sense: 
technological  limitations,  primarily  the  inability  to  avoid  atmospheric  degradation,  will 
require  an  increase  in  laser  output  power  over  current  maximum  operational  levels  in 
order  for  the  weapons  systems  to  be  effectively  lethal. 


It  is  clear  that  tiie  ability  to  assess  the  distribution  of  atmospheric  particulates  on  a 
real-time  basis  is  essential  in  order  to  optimize  the  performance  of  each  electro-optic 
emitter.  Liquid  water  in  the  form  of  cloud  droplets  severely  attenuates  electro-optic 
energy.  In  stratoform  clouds  confined  to  the  lower  atmosphere,  cloud  droplet  size  dis¬ 
tribution  is  a  function  of  the  resident  aerosol  population  and  may  indicate  the  concen¬ 
tration  and  general  composition  of  coincident  aerosols.  Thus,  the  cloud  droplet  size 
distribution  may  serve  as  a  measure  of  the  potential  effectiveness  of  electro-optic  sys¬ 
tems. 

D.  CLIMATIC  VARIABILITY  DUE  TO  THE  RADIATIVE  CONTRIBUTION  OF 

CLOUDS  AND  AEROSOLS 

The  scattering  and  absorption  of  solar  radiation  by  clouds  and  aerosol  particles  in¬ 
fluences  the  earth's  energy  budget  and  thus  can  have  an  impact  on  variations  in  climate, 
fiie  possible  climatic  effects  from  an  increasing  concentration  of  man-made  aerosols 
through  industrial  pollution  and  other  sources  has  received  global  attention  from  the 
scientific  community  and  the  media.  A  description  of  the  effect  of  aerosols  on  the 
transfer  of  radiant  energy  in  the  atmosphere  is  provided  by  Coaklcy  et  al.  (  1983). 
Inclusion  of  the  exact  aerosol  elfect  in  climate  models  is  difficult  since  the  shape,  index 
of  refraction,  size  distribution,  concentration  and  spatial  variability  of  the  aerosols  must 
be  known  in  principle,  and  these  parameters  are  either  too  complex  to  describe  or  vir¬ 
tually  impossible  to  quantify  due  to  unknown  init;ul  conditions  and  temporal  variations. 
Most  climate  models  rely  on  broad  approximations  to  describe  the  distribution  of  these 
constituents. 

Apart  from  the  direct  effect  of  aerosols  on  the  transfer  of  solar  radiation,  Twomcy 
ft  al.  (  1984)  showed  that  aerosols  can  affect  the  radiative  properties  of  clouds  and  create 
an  additional  influence  on  the  radiation  balance.  According  to  Coakley  et  al.  ( 1987),  the 
effect  of  aerosols  on  the  earth's  radiation  budget  through  their  influence  on  clouds  may 
be  several  limes  that  of  the  direct  interaction  of  the  aerosol  with  solar  radiation.  Aer¬ 
osols  are  the  primary  source  of  nuclei  around  which  cloud  droplets  and  ice  crystals  form. 
I  hose  cloud  condensation  nuclei  ( <7(7 N),  which  may  be  less  than  o.oium  in  diameter,  can 
rapidly  form  cloud  droplets  of  several  micrometers  in  radius  upon  entering  a  suitable 
supersaturated  environment.  While  the  minute  CCN  have  only  a  weak  elfect  on  the 
transfer  of  radiation,  the  cloud  droplets  which  form  on  them  have  a  major  impact  on  the 
radiation  balance.  Charlson  et  al.  (  IdST)  showed  the  importance  of  phytoplankton- 
produced  dimethylsulpiiide  (DM8)  gas  in  the  marine  environment.  I  he  sulfate  particles 
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that  form  when  DMS  oxidizes  may  be  the  principal  source  of  CCN  in  marine  air.  This 
implies  that  oceanic  biological  variability  may  play  a  major  role  in  marine  cloud  albedo 
changes  and  in  planetary  climatic  trends.  The  inclusion  of  variable  aerosol  contributions 
to  the  radiative  properties  of  cloud  fields,  such  as  continental  pollution  sources  or  in¬ 
creases  in  DMS  concentration,  has  yet  to  be  completely  examined  on  a  regional  or 
global  scale. 

C.  THE  FIRST  ISCCP  REGIONAL  EXPERIMENT  (FIRE) 

FIRE  is  a  multidisciplinary  research  program  in  support  of  the  International  Satel¬ 
lite  Cloud  Climatology  Program  (ISCCP)  to  study  the  roles  that  clouds  play  in  the 
global  climate.  The  objectives  and  milestones  for  FIRE  are  enumerated  by  Randall  et 
al.  (1984)  and  detailed  by  Cox  et  al.  (1987).  The  43  members  of  the  initial  FIRE  science 
experiment  team  conceived  of  the  following  six  general  objectives  for  the  experiment: 

1.  characterization  of  clouds, 

2.  validation  and  improvement  of  satellite  cloud-retrieval  techniques, 

3.  improvement  of  cloud  radiation  models, 

4.  description  of  cloud  space, 'time  statistics, 

5.  improvement  of  cloud  dynamics  models,  and 

6.  validation  and  improvement  of  general  circulation  model  cloud  parameterizations. 

The  data  for  this  thesis  were  collected  during  the  first  marine  stratocumulus  inten¬ 
sive  field  observations  (IFO),  conducted  from  29  June  to  20  July  1987  off  the  southern 
California  coast  as  shown  in  Fig.  1.  This  IFO  was  the  first  of  two  planned  for  marine 
stratocumulus  cloud  field  measurements.  It  included  a  variety  of  meteorological  and 
surface  radiation  measurements  conducted  from  San  Nicolas  Island,  including  an  in¬ 
strumented  19-mcter  tower  for  moisture  and  heat  llux  measurements,  two  tethered  bal¬ 
loons  which  probed  the  marine  atmospheric  boundary  layer  (MABI.)  structure,  a 
microwave  radiometer,  a  Doppler  radar  and  sodar.  and  rawinsondc  releases  from  the  is¬ 
land  itself  and  from  the  research  vessel  POINT  SL  R  operating  upwind  from  the  island. 
Five  aircraft  participated  in  the  measurement  clfort:  the  National  Aeronautics  and 
Space  Administration  (NASA)  high-altitude  ER-2,  the  University  of  Washington  C-131. 
the  National  Center  for  Atmospheric  Research  (NCAR)  Electra,  a  C-130  operated  by 
the  British  Meteorological  Office,  and  a  Piper  Navahoe  provided  by  the  Naval  Ocean 
Systems  Center  (NOSC).  The  individual  capabilities  of  these  aircraft,  as  well  as  the 


S$ 


I 


•s^.v:.^ 


»wv 


4 

iJ 


118*  W 


iao  w 


Fig.  I.  FIRE  IFO  1987  Operations  Area:  Outlined  area  indicates  region  encom¬ 
passing  case  study  satellite  subscenes;  line  segments  are  applicable  case 
study  aircraft  mission  tracks. 
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ocher  ground  and  satellite-based  measurement  systems,  are  described  by  Cox  et  al. 
(1987).  Aircraft  operations  were  coordinated  with  various  satellite  overflights. 

The  satellite  data  utilized  in  this  thesis  are  from  the  National  Oceanic  and  Atmo¬ 
spheric  Administration  (NOAA)  polar  orbiting  platforms'  (NOAA-9  10)  Advanced  Very 
High  Resolution  Radiometer  (AVHRR)  collected  at  the  Scripps'  Satellite  Oceanography 
Facility.  Cloud  microphysical  data  illustrated  herein  were  obtained  by  sensors  aboard 
the  NCAR  Clectra  during  MABL  cloud-penetrating  missions  in  support  of  the  IFO. 
Data  for  use  in  trajectory  analyses  are  provided  by  the  Fleet  Numerical  Oceanography 
Center  (FNOC)  Navy  Operational  Regional  Analysis  and  Prediction  System  (NORAPS) 
forecast  model.  Synoptic  summaries  are  derived  from  National  Meteorological  Center 
(NMC)  regional  optimal  interpolation  (01)  analyses. 

Kloesel  et  al.  (1988)  provides  a  summary  of  FIRE  IFO  1987  operations  and  synoptic 
conditions.  The  data  in  this  summary  are  utilized  in  this  thesis  to  describe  synoptic  ef¬ 
fects  and  to  fix  the  geographic  positions  of  aircraft  flight  tracks. 

D.  THESIS  OBJECTIVE 

The  objective  of  this  thesis  is  to  demonstrate  the  ability  to  deduce  cloud  micro¬ 
physical  parameters  (primarily  cloud  droplet  size  distribution)  of  marine  stratocumulus 
clouds  utilizing  AVHRR  high  resolution  picture  transmission  (1IRPT)  satellite  data. 
The  boundary7  between  air  masses  of  continental  and  marine  origin  are  highlighted; 
however,  additional  phenomena  previously  observed  by  Coakley  et  al.  (1987),  Coakley 
and  Davies  (1986),  and  Rogers  (1987)  are  observed,  such  as  cloud  edge  effects  and  tracks 
resulting  from  the  discharge  of  ship  stack  diluents.  A  multispectral  image  processing 
technique  is  used  which  combines  AVHRR  channels  1,  3  and  4.  This  technique  relics 
on  the  basic  physical  principles  of  radiative  transfer  in  the  atmosphere  rather  than 
complex  image  enhancement  schemes.  Channel  1  provides  cloud  solar  reflectance  in 
visible  wavelengths  which  is  corrected  for  anisotropy  during  processing.  The  thermal 
emission  contribution  to  the  near-infrared  (N1R)  channel  3  is  removed  using  the  thermal 
emission  measured  in  channel  4,  and  the  remaining  rcllectance  value  is  corrected  for 
anisotropy.  In  addition  to  providing  thermal  emission  values  for  channel  3  correction, 
channel  4  temperature  is  utilized  as  a  threshold  to  mask  the  contribution  from  any  cloud 
elements  which  might  exist  above  the  MAUL.  The  resulting  rcllectance  values  are 
compared  to  theoretical  calculations  and  aircraft  measurements  in  order  to  determine  the 
relationship  between  solar  rcllectance  and  the  microphysical  properties  of  stratocumulus 
clouds  and  their  co-existent  MABL  during  tire  IFO.  This  will  be  accomplished  using 


relationships  between  reflectance  and  cloud  microphysics  illustrated  in  this  thesis.  The 
ability  to  discern  boundary  layer  air  mass  characteristics  through  remote  sensing  tech¬ 
niques  provides  the  opportunity  to  determine  broad-scale  effects  of.MABL  composition 
changes  on  the  earth's  radiative  balance,  as  well  as  assess  the  range-dependent  per¬ 
formance  of  directed  energy  systems. 

In  Chapter  II  the  necessary  theoretical  background  is  presented.  A  discussion  of  the 
relationship  between  the  distribution  of  aerosols,  cloud  droplets  and  droplet  spectra  ob¬ 
served  in  marine  stratocumulus  clouds  is  followed  by  an  illustration  of  applicable 
radiative  transfer  principles,  including  variations  expected  in  marine  stratocumulus 
clouds.  Model-generated  reflectances  for  marine  stratocumulus  clouds  are  illustrated 
and  discussed.  Chapter  HI  describes  the  satellite  data  and  details  image  processing. 
Finally,  Chapter  IV  illustrates  four  case  studies  from  FIRE  IFO  1987.  Each  case  study 
includes  a  synoptic  description,  applicable  satellite  images,  an  analysis  of  cross-sections 
of  solar  reflectance,  a  comparison  with  theoretical  expectations,  and  a  comparison  with 
the  NCAR  Electra  cloud  microphysics  data  for  the  satellite  passes  with  aircraft  mission 
coincidence. 


II.  THEORETICAL  BACKGROUND 


A.  CLOUD  MICROPHYSICS 

This  section  describes  basic  cloud  microphysical  parameters  which  drive  variations 
in  radiative  properties.  Typical  aerosol  distributions  are  examined,  including  their  role 
as  CCN  and  their  variation  between  air  masses  of  continental  and  marine  origin.  This 
is  followed  by  a  discussion  of  the  dependence  of  cloud  droplet  distribution  on  the  initial 
aerosol  distribution.  Finally,  typical  droplet  spectra  in  marine  stratocumulus  clouds  are 
observed. 

1.  Aerosols  and  Aerosol  Distribution  Spectra 

The  atmospheric  particles  that  are  larger  than  the  molecular  scale  can  be  sepa¬ 
rated  into  two  categories:  (1)  cloud  droplets  and  ice  crystals,  and  (2)  aerosols.  Aerosols 
range  in  size  from  about  0.00 l^m  for  small  ions  made  up  of  charged  molecular  clusters, 
to  more  than  10^m  for  the  largest  dust,  salt  and  combustion  particles.  The  radii  of 
aerosols  of  optical  importance  are  between  0.01  and  lChnn;  cloud  droplets  are  much 
larger  and  vary  between  1  and  100^m.  Aerosols  are  concentrated  in  the  lowest  few 
kilometers  of  the  troposphere.  According  to  Chahine  et  ai  (1983)  their  small  mass  and 
size  allow  them  to  remain  suspended  in  the  atmosphere  for  lifetimes  on  the  order  of  a 
week.  Variations  in  concentration  and  dispersal  are  governed  by  wind  patterns,  so  that 
during  their  lifetime  aerosols  may  be  advccted  thousands  of  kilometers  and  chemically 
and/or  physically  interact,  resulting  in  complex  shapes  and  chemical  compositions. 
Brock  (1972)  states  that  nearly  75%  of  the  aerosol  mass  in  the  atmosphere  can  be  at¬ 
tributed  to  natural  and  anthropogenic  primary  sources  such  as  sea  spray  (40%),  wind- 
generated  dust  (20%),  fires  (10%),  and  industrial  operations  and  combustion  products 
(5%).  The  remaining  25%  is  accounted  for  by  photochemical  and  other  chemical 
processes  which  involve  conversion  of  atmospheric  gaseous  constituents  such  as  DMS 
to  small  particles.  Aerosols  continually  undergo  a  variety  of  chemical  and  physical 
transformations,  including  coagulation,  scavenging,  precipitation  washout, 
condensation,  sedimentation,  dispersion  and  mixing.  Slinn  (1975)  has  provided  an  ide¬ 
alized  presentation  of  typical  aerosol  spectra,  presented  in  Fig.  2. 

2.  Effect  of  Aerosols  on  Cloud  Droplet  Formation  and  Distribution 

In  general,  aerosols  can  be  classified  as  hygroscopic,  neutral,  or  hydrophobic 
depending  on  their  affinity  for  water.  This  affinity  depends  on  aerosol  solubility  which 
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Fig.  2.  Idealized  Aerosol  Distribution  Spectra:  Includes  typical  ground-level 
background  distributions  and  the  general  dependence  of  the  spectra  on 
height,  wind  speed,  distance  from  source,  and  surface  heating  (from  Slinn, 
1975). 

allows  ions  to  collect  on  the  surface  of  the  growing  droplet,  attracting  additional  water 
vapor  and  adding  to  the  droplet  mass.  Cloud  droplets  form  most  efficiently  on  the 
hygroscopic  nuclei.  Cloud  droplet  nucleation  on  a  neutral  aerosol  requires  about  the 
same  level  of  supersaturation  as  spontaneous  or  homogeneous  nucleation  (from  Rogers, 
1979;  the  supersaturation  required  is  approximately  300-400  percent  excess  relative  hu¬ 
midity,  and  these  levels  are  never  observed  in  the  atmosphere).  Hydrophobic  aerosols 
require  even  higher  supersaturation  levels.  Soluble  hygroscopic  aerosols  require  much 
lower  levels  of  supersaturation,  although  the  small  aerosols  made  up  of  ion  clusters 
provide  little  improvement  over  homogeneous  nucleation.  and  large  ( lOfim)  aerosols  are 
of  limited  importance  because  of  their  short  residence  time  and  their  relatively  low  con¬ 
centration.  In  the  atmosphere,  the  only  condensation  nuclei  (CN)  that  enter  into  the 
cloud  droplet  formation  process  nucleate  at  supersaturations  of  approximately  1%  or 
less.  This  subset  of  the  CM  distribution  defines  the  CCN.  CN  and  CCN  are  generally 
classified  according  to  their  geometric  size  as  Aitken  or  small  (radius  (r)  <  O.Imiti),  large 


8 


(0.  l^m  <  r  <  l^m)  and  giant  (r  >  l^m).  Aitken  nuclei  are  products  of  combustion  and 
natural  processes  in  the  atmosphere,  while  large  and  giant  nuclei  are  mostly  salt  particles 
that  result  from  the  bursting  of  ocean  surf  bubbles. 

It  is  the  role  of  aerosols  as  CCN  that  allows  the  inference  of  aerosol  distribution 
based  upon  the  resident  cloud  droplet  distribution.  Twomey  and  Warner  (1967)  made 
direct  measurements  of  cloud  droplets  and  nuclei  in  non-precipitating  cumuli  with  tops 
warmer  than  0°C.  They  found  that  cloud  nuclei  concentration  can  be  used  to  infer  the 
cloud  droplet  concentration.  The  droplet  concentration  of  maritime  clouds  corresponds 
to  the  number  of  large  salt  nuclei  observed  when  wind  conditions  are  moderate,  on  the 
order  of  a  few  tens  per  cm3  .  Conversely,  clouds  which  form  under  the  influence  of 
continental  air  usually  contain  several  hundred  drops  per  cm3  while  the  number  of  large 
nuclei  composed  of  sea  salt  remain  at  about  10  per  cm 3  .  The  continental  aerosol  dis¬ 
tribution  also  contains  significant  quantities  of  soluble  sulfate  components  such  as 
sulfuric  acid  or  ammonium  sulfate  from  various  natural  and  man-made  sources,  and  in¬ 
soluble  particles  from  the  soil.  Charlson  ec  al.  (2987)  indicates  that  the  droplet  concen¬ 
tration  of  maritime  clouds  may  be  more  closely  tied  to  the  ambient  source  level  of 
D MS-induced  sulfates.  Low  droplet  concentrations  in  marine  clouds  may  result  from  a 
lack  of  suitable  CCN.  The  relative  importance  of  each  of  the  aerosol  constituants  de¬ 
pends  upon  the  spatial  and  temporal  history  of  the  associated  air  mass  and  its  inter¬ 
action  with  land  and  sea  surfaces. 

An  example  of  a  continental  influence  on  aerosol  distribution  and  its  associated 
effect  on  the  concentration  of  cloud  droplets  is  illustrated  by  Warner  and  Twomey 
(1967).  An  additional  source  of  CCN  observed  over  maritime  regions  east  of  Australia 
was  isolated  as  cane  lire  smoke  products  resulting  from  inland  burning  during  the  har¬ 
vest  season.  In  this  CCN-rich  air,  cloud  droplet  concentrations  in  low-level,  warmer 
than  0°C  clouds  were  observed  in  concentrations  5  to  10  times  higher  than  had  previ¬ 
ously  been  observed  in  the  same  maritime  region,  and  10  or  more  times  higher  than  the 
median  concentration  expected  in  maritime  cumuli.  Hindman  et  al.  (1977)  observes  the 
dependence  of  cloud  droplet  size  distribution  on  size-catcgorizcd  CCN  in  warm,  non¬ 
raining  clouds.  The  results  of  both  measurements  and  calculations  show  that  high  con¬ 
centrations  of  large  and  giant  CCN  lead  to  droplet  size  distributions  which  contain  large 
droplets  only  if  low  concentrations  of  small  CCN  arc  present.  If  high  concentrations 
of  small  CCN  are  present,  the  droplet  size  distributions  contain  few,  if  any,  large  drop¬ 
lets,  regardless  of  the  concentration  of  the  larger  CCN.  High  droplet  concentrations  are 


commonly  observed  in  clouds  which  contain  few  large  droplets,  while  low  droplet  con¬ 
centrations  usually  indicate  the  presence  of  a  greater  percentage  of  large  droplets.  It 
appears  that  the  relative  concentration  of  small  and  large  CCN,  not  the  concentration 
of  giant  CCN,  determines  the  extent  of  large  droplets  in  low-level  warm  cumuli. 

3.  Cloud  Droplet  Spectra  Observed  in  Marine  Stratocumulus  Clouds 

The  droplet  spectra  in  marine  stratocumulus  clouds  is  largely  dependent  upon 
the  number  and  chemical  nature  of  the  CCN,  which  are  in  turn  controlled  by  air  mass 
source,  sea-sprav  aerosols  as  regulated  by  the  surface  wind  speed,  and  oceanic  biological 
activity.  Noonkester  (1984)  made  detailed  observations  of  cloud  droplet  spectra  in  ma¬ 
rine  stratus  layers  110-150  kilometers  southwest  of  San  Diego,  CA,  in  the  spring  anu 
summer  of  1981  (this  location  coincides  with  the  FIRE  IFO  1987  operating  area).  Two 
days  were  selected  for  study:  one  day  reflecting  the  influence  of  a  marine  air  mass  (29 
May  1981)  and  the  other  representing  the  influence  of  continental  air  (18  August  1981). 
These  associations  were  based  on  estimates  of  air  parcel  trajectories  from  synoptic  ana¬ 
lyses.  It  was  observed  that  droplet  spectra  on  the  continental  and  marine  days  were 
different,  primarily  due  to  the  existence  of  many  more  small  droplets  in  air  of  continental 
origin.  In  the  upper  third  of  cloud  layers,  the  continental  regime  surrendered  values  of 
500-600  droplets  per  cm 3  at  an  average  modal  droplet  radius  (i.e.,  the  radius  of  peak 
concentration)  near  4^m.  The  marine  case  gave  droplet  concentrations  of  200-300  per 
cm 3  at  an  average  modal  radius  near  6/un.  Fig.  3  shows  measurements  of  these  param¬ 
eters  and  the  mean  droplet  radius  as  a  function  of  the  distance  from  the  cloud  base  for 
both  air  mass  types.  Values  at  various  flight  levels  within  the  cloud  are  listed  in 
Table  1. 

The  equivalent  liquid  water  content  (EQLWC)  is  defined  as  the  mass  of  the 
liquid  water  contained  in  all  of  the  cloud  droplets  or  ice  crystals  per  unit  volume.  The 
EQLWC  measured  in  the  vertical  by  Noonkester  (1984)  reveal  an  average  value  of  0. 1 
g/m3  at  mid-cloud  levels  increasing  to  0.3  g/m3  near  the  cloud  top  for  both  continental 
and  marine  cases  as  shown  in  Fig.  4. 

Noonkester  (1984)  also  observed  that  the  aerosol  spectra  in  air  of  continental 
origin  may  vary  based  upon  the  length  of  time  air  parcels  may  be  modified  by  the  marine 
environment  as  they  pass  over  ocean  surfaces.  An  air  mass  off  the  southern  California 
coast  may  consist  mainly  of  "fresh"  continental  CCN  due  to  offshore  flow  over  the  ad¬ 
jacent  coastline,  or  may  contain  "aged"  CCN  advected  into  the  region  from 
coastal, continental  sources  farther  north  (the  Pacific  Northwest  or  possibly  Canada). 
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Fig.  3.  Profiles  of  Stratocumulus  Cloud  Droplet  Parameters:  Samples  are  from 
continental  air  (18  Aug)  and  marine  air  (29  May)  through  the  cloud  layer. 
Dots  are  sampling  elevations  (from  Noonkester,  1984). 

The  CCN  in  the  latter  have  been  modified  through  the  addition  of  marine  aerosols  and 
the  loss  processes  described  in  the  previous  section. 

B.  RADIATIVE  TRANSFER  PRINCIPLES 

The  satellite  data  utilized  in  this  thesis  comprise  three  of  the  five  AVHRR  channels 
employed  on  the  NOAA-9/10  polar  orbiting  platforms.  These  include  channel  1 
(0.58-0.68/im  bandwidth;  centered  at  0.63Mni)  in  visible  wavelengths,  channel  3 
(3. 55-3. 93, um;  3.70Mm)  in  the  NIR,  and  channel  4  (10.30-1  lJO^m;  11.00/im)  in  the  in¬ 
frared  (IR)  region.  The  following  section  provides  a  discussion  of  the  radiative  transfer 
characteristics  in  these  portions  of  the  electromagnetic  spectrum  in  order  to  understand 
the  methods  developed  and  utilized  for  satellite  image  processing.  Where  appropriate, 
simplifications  are  made  based  on  the  assumption  that  only  satellite  pixels  encompassing 
stratocumulus  clouds  will  be  considered  in  the  radiative  balance. 

1.  Total  Radiance 

The  total  satellite-measured  spectral  radiance  in  units  of  W/m^-sr-zirn  is  given 
by 
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Fig.  4. 


Profiles  of  Equivalent  Liquid  Water  Content  (EQLWC):  The  solid  curve 
represents  the  marine  regime,  while  the  dashed  curve  represents  the  conti¬ 
nental  regime.  The  zero  reference  level  is  cloud  base,  defined  to  be  where 


EQLWC  is  0.02  g!mi  (from  Noonkester,  1984). 


TABLE  1 

CONTINENTAL/MARINO  STRATOCL'ML'LL'S  CLOL'D  PARAMETERS 

(alter  Noonkester,  19S4) 

Number  of 

Mean 

z' 

EQLWC 

Droplets 

Droplet 

,  lx 

at  r 

Radius  r 

Level 

(m  from 

{glut) 

cloud  base) 

(cm~3) 

(M"  t) 

Marine  Air  M ass  Regime 

1 

-380 

0.00103 

21 

0.87 

4 

-63 

0.00314 

194 

0.83 

8 

18 

0.0337 

260 

2.00 

11 

167 

0.217 

296 

4.82 

12  (CLD  TOP) 

274 

0.288 

22S 

13 

301 

o.ooool 

2 

Continental  Air  Mass  Regime 

1 

-407 

0.00014 

in 

5 

-62 

0.00054 

8 

17 

0.0374 

12 

121 

0.199 

14  (CLD  TOP) 

222 

0.263 

575 

15 

254 

0.00002 

90 

0.31 

The  following  definitions  apply: 

L, 

= 

upwclling  radiance, 

/ 

= 

wavelength. 

0 

= 

satellite  zenith  angle. 

<t> 

= 

azimuth  angle, 

U 

= 

surface  radiance. 

6 

= 

optical  depth. 

- 

= 

vertical  distance. 

- 

beam  attenuation  cocllicient. 

t 

= 

emissivity, 

X 

= 

elemental  position  vector. 

D 

= 

Planck  function, 

<•% 

= 

volume  scattering  coellicicnt. 

/*(»// s)  =  scattering  phase  function, 

L  =  incident  radiance, 

d,  =  scattering  optical  depth,  and 

(?„  =  solar  zenith  angle. 

The  first  term  on  the  right  hand  side  of  Eq.  1  represents  the  direct  surface 
transmittance,  that  is,  the  probability  of  surface  reflected  or  emitted  radiant  energy 
passing  unobstructed  through  the  atmosphere.  The  contribution  to  total  radiance  from 
this  term  is  small  for  satellite  pixels  filled  with  clouds  of  sufficient  thickness.  The  surface 
contribution  is  masked  by  the  clouds  in  visible  wavelengths  and  is  dominated  by  the 
clouds  in  the  IR.  The  reflection  and  emission  from  cloud  surfaces  in  the  NIR  comprise 
virtually  all  of  the  upwelling  radiation  in  these  pixels.  In  addition,  reflection  off  the 
ocean  surface  is  small  in  the  visible  and  NIR  away  from  regions  of  sunglint.  In  the  case 
of  IR  wavelengths,  the  desired  quantity  is  the  ncar-blackbody  emission  of  the 
stratocumulus  clouds.  Since  the  emissivity  (absorptivity)  of  these  clouds  is  close  to  1.0, 
the  clouds  effectively  cancel  any  contribution  from  below  the  cloud  layer. 

The  second  term  on  the  right  hand  side  of  Eq.  1  accounts  for  radiance  sources 
and  sinks  in  the  path  between  the  surface  and  the  satellite.  If  it  is  assumed  that  the 
transmissivity  between  the  stratocumulus  cloud  top  and  the  satellite  is  1.0  (no  additional 
clouds  or  scattering/absorbing  media  above  the  stratocumulus  cloud  top),  then  this  term 
includes  radiant  energy  scattered  out  of,  and  into,  the  sensed  path  by  elements  within  the 
cloud,  absorption  in  the  cloud  which  removes  radiation  from  the  sensed  path,  and  cloud 
droplet  emission  into  the  path.  The  assumption  of  perfect  transmissivity  above  the 
clouds  is  valid  since  the  attenuation  contribution  of  aerosols  above  the  MABL  is  small. 
Also,  water  vapor  content  falls  off  sharply  above  the  MABL,  especially  in  the  cases 
during  FIRE  IFO  19S7  where  the  region  remained  under  the  influence  of  strong,  drying 
subsidence.  This,  in  conjunction  with  the  fact  that  channels  1  and  3  are  in  atmospheric 
"window"  regions,  minimizes  the  effect  of  water  vapor  on  transmissivity.  These  same 
assumptions  arc  deemed  valid  by  Bell  and  Wong  (1981)  when  describing  the  NIR  radi¬ 
ation  received  by  satellite  sensors. 

With  the  above  assumptions  in  effect,  Eq.  1  may  be  simplified  to  the  form 

/.,(/.  0,  <(>)  =  £(/)/?(/, T)  +  yr( / .  00,  0,  «/i )L  cos  00,  (5) 

where  >■,(/,  0,,,0,<t>)  is  the  directional  reflectance  of  the  cloud  surface  received  by  the 
satellite.  The  first  term  on  the  right  hand  side  of  Eq.  5  is  the  contribution  to  measured 


radiance  due  to  thermal  emission  as  determined  by  the  Planck  function  (D)  and  the 
emissivitv  (c)  of  the  cloud  surface.  The  Planck  function  relates  emitted  monochromatic 
intensity  with  wavelength  and  cloud  surface  absolute  temperature  as  follows: 


2  he' 


j  5^(hcjJ.kl)  _  j 


(6) 


with 

h  =  Planck's  constant  (6.6262(10-34)  J  sec), 

c  =  velocity  of  light  (2.998(10*)  m,  sec), 

k  =  Boltzmann's  constant  ( 1.3806(ICH3)  J,  deg),  and 

T  =  brightness  temperature. 

The  second  term  on  the  right  hand  side  of  Cq.  5  is  the  contribution  to  the  up- 
welling  radiance  due  to  solar  reflection  as  determined  by  the  incident  solar  radiation 
(L)  weighted  by  the  cosine  of  the  solar  zenith  angle  ((>„)  and  the  directional  reflectance 
(•;,).  The  directional  reflectance  describes  reflected  energy  from  some  incident  direction 
into  a  particular  outward  direction  as  defined  by  by  the  solar  zenith  angle,  the  satellite 
zenith  angle  (0),  and  the  relative  azimuth  (</>).  Relative  azimuth  is  defined  as  the  angle 
between  the  projection  of  the  sun-to-pixel  and  satellite-to-pixel  line  segments  on  a  hori¬ 
zontal  plane. 

2.  Cloud  Radiative  Processes 
a.  Overview 

Solar  radiation  which  encounters  atmospheric  particulates  (molecular 
constituants,  aerosols,  and  cloud  droplets  or  ice  crystals)  may  be  transmitted  through  the 
atmosphere,  or  may  be  absorbed  or  scattered  by  those  particulates.  If  the  scattering 
particles  arc  much  smaller  than  the  wavelength  of  the  radiant  energy,  the  interactions 
are  described  by  Rayleigh  scattering.  Rayleigh  scattering  efficiency  is  inversely  propor¬ 
tional  to  the  fourtii  power  of  wavelength  and  as  such  is  not  important  for  the  wave¬ 
lengths  and  droplet  sizes  of  interest  in  this  thesis.  If  the  particle  radius  and  wavelength 
are  comparable,  the  scattering  process  is  described  by  Mie  theory  (I. ion,  1980).  I'nder 
Mie  theory,  the  volume  scattering  coefficient  (rrj,  in  units  of  inverse  length,  describes  the 
total  scattering  in  a  thin  horizontal  section  of  the  atmosphere  and  is  defined  as  follows: 


Fig.  5.  Mie  scattering  efficiency  factor  (Q):  Illustrated  as  a  function  of  the  size 
parameter  x  =  Inr/j.  .  The  real  refractive  index  is  1.33,  with  results  shown 
for  four  values  of  the  imaginary  part  of  the  refractive  index  (after  Hansen 
and  Travis,  1974). 

where  m  is  the  index  of  refraction.  The  Mie  scattering  efficiency  factor  (Q)  is  illustrated 
in  Fig.  5  for  variable  indices  of  refraction.  The  real  part  of  the  complex  index  of  re¬ 
fraction  for  pure  water  is  1.33,  while  the  imaginary  part  ranges  from  0  for  a  perfect  re¬ 
flector  to  positive  numbers  less  than  unity  for  absorbing  media.  The  major  maxima  and 
minima  in  the  curves  are  due  to  interference  of  radiant  energy  diffracted  and  transmitted 
by  a  particle,  while  the  rapid  fluctuations  result  from  radiant  paths  which  graze  the  edge 
of  the  particle.  Note  that  the  Mie  scattering  efficiency  asymptotically  approaches  a 
value  of  2  for  relatively  large  particle  sizes.  Thus,  geometric  optics,  which  describe 
scattering  mechanisms  where  the  particles  are  much  larger  than  the  wavelength,  define 
a  cross  sectional  area  approximately  twice  that  of  a  spherical  particle  (from  Liou.  1980) 
within  which  scattering  will  occur.  Using  Eq.  7  to  describe  the  scattering  processes 
within  a  cloud,  it  is  clear  that  the  effects  of  cloud  droplets,  being  much  larger  than  their 
associated  CCN,  dominate  the  radiative  scattering.  The  volume  scattering  coefficient 


may  typically  be  as  much  as  100  times  as  large  for  cloud  droplets  as  for  the  CCN  based 
on  their  comparative  sizes.  This  fact  allows  the  contribution  from  CCN  scattering  to  the 
total  cloud  reflectance  to  be  ignored. 

Reflectance  from  a  cloud  for  the  wavelengths  of' interest  depends  upon  these 
Mie  scattering  mechanisms.  The  scattering  optical  depth  (Js)  provides  a  measure  of  the 
scattering  based  upon  the  physical  depth  of  the  scattering  layer  (e.g.,  in  this  thesis,  the 
stratocumulus  cloud  thickness)  as  follows: 


Cloud  reflectivity  is  also  a  function  of  the  angular  variation  of  scattering  for 
each  single  scattering  event.  This  variation  is  the  scattering  phase  function  (Tf/h))  and 
is  used  to  describe  tendencies  toward  directionally  dependent  scattering.  Mie  scattering 
phase  functions  are  dominant  in  the  forward  direction  (i.e.,  more  radiance  is  scattered 
in  the  direction  of  the  original  ray  than  in  the  reverse  direction),  while  Rayleigh  scatter¬ 
ing  is  nearly  the  same  forward  and  back.  The  properties  vary  with  wavelength  and 
particle,' droplet  size.  In  the  case  of  Mie  scattering,  as  the  particles  or  droplets  increase 
in  size,  more  forward  scattering  occurs.  Hunt  (1972)  illustrates  phase  function  diagrams 
for  size  distributions  of  cloud  droplets  and  ice  crystals  found  in  various  cloud  types,  and 
confirms  that  a  decrease  in  cloud  droplet  modal  radius  results  in  greater  back-scattering 
of  radiation  by  the  droplets.  While  the  reflectance  from  a  cloud  is  dependent  on  the 
nature  of  single  scattering  characteristics,  it  must  be  remembered  that  radiant  energy 
from  a  cloud  source  at  visible  wavelengths  will  have  been  scattered  many  times,  effec¬ 
tively  smoothing  out  the  complex  angular  patterns  observed  in  single  scattering  events. 
In  the  NIR,  the  probability  of  absorption  limits  the  number  of  individual  scattering 
events  so  that  more  of  the  directional  scattering  defined  by  the  phase  function  is  main¬ 
tained. 

The  ratio  of  the  volume  scattering  coefficient  to  the  beam  attenuation  co¬ 
efficient  (ojt 7e)  is  called  the  single  scattering  albedo  (<ou)  and  describes  the  probability 
of  scatter  for  a  single  interaction  between  radiant  energy  and  a  particle  or  droplet.  The 
single  scattering  albedo  varies  with  wavelength  and  particle,  droplet  composition. 

As  illustrated  by  Hunt  (1972).  changes  in  cloud  radiative  properties  for  the 
wavelengths  of  interest  arc  consistent  with  changes  in  cloud  droplet  size  distribution  and 
concentration,  and  scattering  optical  depth  which  is  a  function  of  cloud  thickness.  It 


was  shown  in  a  previous  section  that  the  cioud  droplet  distribution  depends  upon  the 
initial  distribution  of  CCS  present  in  the  air  mass.  It  follows  logically  that  changes  in 
the  reflective  properties  of  clouds  may  be  due  to  variations  in  aerosol,  CCN  distribution. 
This  relationship  is  illustrated  by  Twomey  and  Cocks  (19S2)  and  discussed  by  Coakley 
el  al.  (1987)  and  Wetzel  and  Vonder  I  laar  ( 1986).  Twomey  ( 1977)  observed  that  in  most 
clouds  the  injection  of  an  additional  source  of  small  CCN  will  increase  the  cloud 
reflectivity.  Hansen  and  Pollack  (1970)  found  that  the  volume  extinction  coefficient,  the 
single  scattering  albedo,  the  shape  of  the  phase  function,  and  the  asymmetry  factor  (i.e., 
the  first  moment  of  the  phase  function;  in  Mie  scattering  it  denotes  the  relative  strength 
of  forward  scattering)  all  depend  mainly  on  cloud  droplet  modal  radius  and  not  on  var¬ 
iations  in  the  shape  of  the  size  distribution.  Cloud  reflectance  will  also  increase  if  its 
scattering  optical  depth  increases  without  a  compensating  gain  in  absorptivity.  Such  an 
effect  occurs  when  there  is  an  increase  in  EQLWC,  although  the  effect  varies  with 
wavelength  and  cloud  thickness  as  will  be  discussed  in  the  next  section. 
b.  Spectral  Variations 

Channel  1  radiance  is  composed  of  reflected  solar  radiation.  Stratocumulus 
clouds  in  this  channel  will  reflect  10  to  15  times  as  much  radiation  than  the  background 
ocean  surface  (discounting  sunglint).  Channel  3  radiance  is  composed  of  a  combination 
of  reflected  and  emitted  radiant  energy.  According  to  Coakley  et  al.  (1987)  the  scattering 
cross  section  of  a  cloud  droplet  at  this  wavelength  is  approximately  proportional  to  its 
geometric  cross  section  while  the  absorption  cross  section  is  approximately  proportional 
to  its  volume.  Thus  the  ratio  of  scattering  to  absorpt  ion  increases  as  droplet  radius  de¬ 
ceases.  The  size  dependence  on  absorption  allows  the  volume  absorption  coefficient 
(aj  and  d’-op'et  radius  to  be  related  to  the  single  scattering  albedo  (from  Twomey  and 
Ox\  -v 2 )  by 

l  -  w0  oe  <7ar .  (9) 

The  radiative  properties  of  stratocumulus  clouds  in  channel  3  are  sensitive  to  changes  in 
cloud  droplet  radius.  Coakley  and  Davies  (19S6)  indicates  that  the  sensitivity  of 
reflectance  to  droplet  size  is  greatest  in  channel  3  due  to  the  absorption  occurring  at  this 
wavelength.  The  additional  scattering  in  conjunction  with  a  decrease  in  droplet  radii 
leads  to  increased  cloud  reflectivity.  Thus,  an  injection  of  CCN -rich  air  into  a  particular 
cloud  field  and  the  associated  increase  in  droplet  concentration  would  be  compensated 
by  a  shift  toward  smaller  droplets  if  the  EQLWC  remains  unchanged.  Smaller  cloud 


droplet  distributions  show  increased  reflectance  in  channel  3  as  a  result  of  increased 
back-scattering,  greater  Mie  scattering  efficiency,  and  the  change  in  the  ratio  of  scatter¬ 
ing  to  absorption.  Higher  droplet  concentrations  and  increased  back-scattering  result 
in  greater  channel  1  reflectance.  In  cloud  regions  where  the  shift  to  smaller  droplet  size 
distribution  is  accompanied  by  a  loss  of  EQLWC  and/or  optical  depth,  the  increase  in 
channel  3  reflectance  remains  prevalent  while  channel  1  reflectance  may  actually  de¬ 
crease.  Arking  and  Childs  (1985)  demonstrates  the  ability  to  deduce  general  micro¬ 
physical  model  parameters  for  cloud  droplets  from  the  total  channel  3  radiance; 
however,  the  ability  of  this  technique  to  resolve  precise  cloud  droplet  size  distribution 
information  is  limited  by  the  variable  contribution  of  thermal  emission  in  this  channel. 
Hunt  (1972)  and  Hansen  and  Pollack  (1970)  report  values  for  solar  reflectance  for  both 
visible  and  N1R  wavelengths.  These  values  show  that  in  channel  1,  reflectance  is 
strongly  dependent  upon  cloud  thickness  and  to  a  lesser  degree  the  EQLWC,  the  modal 
radius,  and  the  solar  zenith  angle.  In  channel  3  the  dependencies  on  EQLWC  and  cloud 
thickness  are  not  as  great  as  the  dependence  on  solar  zenith  angle  and  modal  radius. 
The  single  scattering  albedo  in  channel  3  for  stratocumulus  cloud  droplets  has  an  aver¬ 
age  value  near  0.89;  therefore,  there  is  a  significant  probability  for  absorption  after  rel¬ 
atively  few  scattering  events.  As  a  result,  these  clouds  in  channel  3  essentially  have  an 
infinite  optical  depth  with  a  cloud  thickness  of  just  a  few  tens  of  meters.  Increases  in 
EQLWC  for  clouds  of  sufficient  thickness  do  not  affect  the  scattering  optical  depth,  so 
there  is  no  associated  increase  in  reflectivity. 

Essentially  all  of  the  radiant  energy  in  channel  4,  and  the  radiance  in  chan¬ 
nel  3  which  is  not  due  to  reflection,  is  emitted  thermal  radiation.  Assuming  that 
stratocumulus  clouds  are  optically  thick  and  reflectance  is  discounted,  the  emissivity  in 
channel  3  can  be  determined  by  using  the  brightness  temperature  obtained  from  channel 
4  measurements  and  the  Planck  function.  The  effects  of  solar  reflection  on  brightness 
temperature  can  be  neglected  in  channel  4  since  the  emissivity  of  stratocumulus  clouds 
of  applicable  optical  depth  is  very  close  to  1.0.  With  the  contribution  from  reflectance 
removed  from  channel  3,  the  remainder  of  the  upwelling  radiance  is  thermal  emission. 
As  will  be  seen  in  Chapter  III,  clouds  other  than  those  which  exist  in  the  MAUL  will 
be  identified  utilizing  channel  4  brightness  temperature  so  that  their  more  complex 
absorption,  reflection  characteristics  are  not  a  factor. 


c.  Anisotropy 

Taylor  ami  Stowe  (1984)  showed  that  hemispheric  reflectance  is  related  to 
the  directional  reflectance  by  the  anisotropic  reflectance  factor  (;r)  by 


P(A)  = 


0(l.  0 ,  f) 


(10) 


The  dependence  of  reflectance  on  solar  zenith  angle,  satellite  zenith  angle,  and  relative 
azimuth  was  observed  by  Ruff  and  Gruber  (1983)  by  noting  the  variation  in  reflectance 
on  an  ocean  surface  at  the  same  geographic  point  over  subsequent  satellite  orbits.  The 
anisotropic  reflectance  factors  were  derived  by  Taylor  and  Stowe  (1984)  from 
N1MBLS-7  Earth  Radiation  Budget  sensor  data  and  represent  spectrally  integrated 
albedo  between  0.2/im.  and  50^m.  They  are  categorized  into  discrete  angular  bins,  each 
of  which  span  approximately  12-15  degrees  of  satellite  zenith  angle  and  10-20  degrees 
of  relative  azimuth.  These  bins  are  separated  into  ten  categories  of  solar  zenith  angle 
with  intervals  of  6  to  25  degrees.  Factors  are  listed  for  eight  different  surfaces:  clear 
land,  snow,  ice,  ocean,  high  ice  clouds,  and  high,  middle,  and  low  water  clouds. 

The  amount  that  the  anisotropic  reflectance  factor  is  different  from  a  value 
of  1.0  gives  the  relative  importance  of  the  factor  in  reflectance  calculations.  For  exam¬ 
ple,  a  factor  of  0.75  indicates  that  the  directional  reflectance  is  less  than  the  isotropic 
reflectance  by  25%.  The  factors  for  low  water  cloud  (utilized  in  this  thesis)  range  from 
about  0.5  to  2.5  for  the  range  of  applicable  solar  zenith  angles  (30-70  degrees).  In  gen¬ 
eral,  low  cloud  surfaces  appear  brighter  when  viewed  from  overhead  than  when  viewed 
from  an  angle  near  the  horizon.  The  adaptation  of  anisotropic  reflectance  to  the  satellite 
data  processing  scheme  will  be  discussed  in  Chapter  III. 

C.  MODEL  REFLECTANCE  OF  MARINE  STRAFOCUMULUS  CLOUDS 
1.  Reflectance  Calculations 

To  arrive  at  theoretical  estimates  of  cloud  reflectance,  realistic  cloud  droplet 
distributions  must  be  generated  which  represent  anticipated  conditions  in  the  marine 
stratocumulus  environment.  Model  cloud  droplet  distributions  arc  generated  using  the 
modified  gamma  distribution  (from  Dcirmendjian,  1969)  as  follows: 


where  the  parameters  (A),  (a),  and  (y)  are  empirically  derived  constants.  Values  of  A 
adjust  the  amplitude  of  each  cloud  droplet  size  distribution  based  upon  specified  modal 
radius  and  EQLWC,  and  do  not  affect  distribution  shape.  Typical  values  for  these 
constants  arc  specified  by  Welch  et  al.  (19S0)  and  a  summary  of  those  utilized  in  this 
thesis  are  listed  in  Table  2. 

Modal  mdii  (rj  of  -him  and  8/rm  were  selected  to  represent  realistic  continental 
and  marine  air  masses,  respectively.  Two  typical  levels  of  EQLWC  were  selected  (0.4 
glm*  and  0.8 g/m-')  to  observe  variations  in  reflectance  based  on  this  parameter.  Values 
for  the  empirical  constants  were  chosen  so  as  to  give  a  typical  droplet  size  distribution 
for  marine  stratocumulus  clouds  near  the  cloud  top.  In  addition,  constants  were  selected 
for  two  other  droplet  size  distributions,  one  strongly  peaked  at  the  modal  radius  and  the 
other  with  a  broad  spread  of  droplet  sizes  in  order  to  observe  the  effect  of  variations  in 
droplet  distribution  on  the  calculated  reflectance.  The  droplet  size  distribution  curves 
for  two  combinations  of  modal  radius  and  EQLWC  are  shown  in  Tig.  6. 


TABLE  2 

MODIFIED  GAMMA  DISTRIBUTION  EMPIRICAL  CONSTANTS 


Using  the  12  resultant  droplet  size  distributions,  radiance  characteristics  for 
each  distribution  were  calculated  utilizing  Mie  code  provided  by  Wiscombe  (1980).  fa¬ 
bles  3  and  4  list  the  results  for  each  case.  The  three  combinations  of  a  and  >■  which  de¬ 
scribe  distribution  shape  are  labeled  as: 


i 

i 


NUMRER/CM**3/M1CR0METER  NUMBER/CM**3/M1CR0METER 

.0  20.0  30.0  40.0  60.0  60.0  0.0  6.0  10.0  16.0  20.0  25.0  30.0 


M0DE-R=8  EQLWC=0.4 


_ RADIUS  (MICROMETERS) _ 

Fig.  6.  Model  Cloud  Droplet  Distributions:  Modal  radius  =  S^m,  EQLWC  = 
0Agjmi  (top)  and  O.Sg/m3  (bottom).  Dl,  D2,  and  D3  are  the  three  modified 
gamma  distributions. 
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TABLE  3 

MIE  OUTPUT  PARAMETERS  FOR  0.63 gun 


rc  =  4.0/a/n 


Distribution  D1 
Distribution  D2 
Distribution  D3 


rc  —  4.0gm 


a  =  2.0  y  =  1.19 

a  =  5.0  y  =  2.41 
a  =  5.0  y  =  1.30 


rc  =  8 .0fx/n 


rc  — 


{m)  EQLWC  —  0.Sgjm'!  EQLIVC  =  OAg/rn  EQLWC  =  O.Sglm1  EQLWC  =  0.4g/m 


D1  D2  D3  D1  D2  D3  D1  D2  D3  D1 


1.0  1.0  1.0 


.857  0.845  0.849 


Scattering  Optical  Depth  (<> ) 


56.0 

37.3 

18.7 

7.5 

3.7 

1.9 

0.7 


Single  Scattering  Albedo  (o>0) 


1.0  1.0  1.0  1.0 


Asymmetry  i7actor  (gf) 


0.845  0.849  0.869  0.861  0.864  0.869  0.861 


85.5 

54.8 

57.0 

36.5 

28.5 

18.3 

11.4 

7.3 

5.7 

3.7 

2.9 

1.8 

1.1 

0.7 

D3 

Dl 

mam 

D3 

) 

0.111 

0.037 

IhxtodBnBBE 

83.3 

55.5 

27.8 

11.1 

5.6 

2.8 

1.1 

_ 

26.3 

15. 3 
9.1 

3.7 

1.8 
0.9 
0.4 

50.5 

33.7 

16.8 

6.7 
3.4 

1.7 
0.7 

41.6 

27.7 
13.9 

5.5 

2.8 

1.4 

0.6 

Dl:  (a  =  2.0,  y=  1.19), 

D2:  (a  =  5.0,  y  =  2.41),  and 
D3:  (a  =  5.0,  y  =  1.30). 

The  delta- Eddington  approximation  lor  multiple  scattering  in  clouds  (from  Joseph  ei  ai, 
1976)  was  applied  at  three  representative  solar  zenith  angles  (30.  45.  and  60  degrees)  and 
cloud  thickness  ranging  from  10  to  750  meters.  The  resulting  reflectance  values  arc  listeJ 
in  Tables  5  and  6. 
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TABLE  4 

MIE  OUTPUT  PARAMETERS  FOR.  3.70/im 


Distribution  D1  =>  a  =  2.0  y  =  1.19 

Distribution  D2  =>  a.  =  5.0  y  =  2.41 

Distribution  D3  =>  oc  =  5.0  y  =  1.30 


(>n) 

rc  =  4.0  jim 

EQLWC  -  Q.Zgltn 

rc  =  4.0  jj.ni 

EQLWC  =  0Aglm'' 

rc  =  S.Oju/n 

EQLWC  =  O.Sg/m1 

rc  =  8.0  nm 

EQLWC  -  OAglm 

D 1 

D2 

D3 

D1 

D2 

D2 

D3 

D 1 

■SB 

1)3 

Volume  Scattering  Coefficient  (<7S) 

0.156 

0.345 

0.259 

0.078 

0. 1 73 

0.135 

0.108 

BUM 

0.05-1 

Scattering  Optical  Depth  (<),) 

750 

117.0 

258.8 

194.3 

58.6 

129.8 

96.8 

50.2 

81.0 

25.0 

50.7 

40.5 

500 

7S.0 

172.5 

129.5 

39.1 

86.5 

64.5 

33.5 

54.0 

16. S 

33.8 

27.0 

250 

39.0 

86.3 

64.8 

19.5 

43.3 

32.3 

16.7 

27.0 

8.4 

16.9 

13.5 

100 

15.6 

34.5 

25.9 

7.8 

17.3 

12.9 

6.7 

13.5 

10.8 

3.4 

6.8 

5.4 

50 

7.8 

17.3 

13.0 

3.9 

8.7 

6.5 

3.4 

6.8 

5.4 

1.7 

3.4 

2.7 

25 

3.9 

8.6 

6.5 

1.9 

4.3 

3.2 

1.7 

3.4 

2.7 

O.S 

1.7 

1.3 

10 

1.6 

3.5 

2.6 

0.8 

1.7 

1.3 

0.7 

1.4 

l.l 

0.3 

0.7 

0.5 

Single  Scattering  Albedo  («>„) 

0.911 

0.954 

0.940 

0.9 1 1 

0.954 

0.940 

0.845 

0.899 

0.883 

0.845 

0.899 

Asymmetry  Factor  (#,) 

0.783 

0.755 

0.750 

0.783 

0.S01 

0.817 

on 

2.  Retlectance  Variability 

Fig.  7  through  9  illustrate  the  dependency  of  reflectance  on  cloud  droplet  size 
distribution.  EQLWC,  and  cloud  thickness  for  the  center  wavelength  of  AVMR.Il  chan¬ 
nels  1  and  3  based  on  delta-Eddington  calculations.  The  expected  increase  in  reflectance 
for  droplet  distributions  of  smaller  modal  radius  is  shown  in  Fig.  7.  The  trends  in 
reflectance  support  the  connection  of  higher  reflectances  in  botii  channels  from  clouds 
under  the  influence  of  higher  CCN  concentrations  when  the  EQLWC  and  cloud  thick¬ 
ness  remains  unchanged.  Fig.  8  reveals  the  increase  in  channel  1  reflectance  with  in¬ 
creasing  EQLWC  while  the  channel  3  rellcctance  remains  constant.  The  change  in 
EQLWC  at  a  constant  modal  radius  changes  the  scattering  optical  depth  of  the  cloud  in 
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channel  1  and  the  related  reflectance  is  modified.  In  contrast,  the  absorption  by  liquid 
water  in  channel  3  causes  the  scattering  optical  depth  to  appear  infinite  in  clouds  of 
appreciable  thickness.  At  cloud  thicknesses  of  approximately  50-100  meters  or  less,  de¬ 
creases  in  EQLWC  and  the  associated  loss  in  optical  depth  begin  to  dominate  micro- 
physical  effects.  The  channel  3  reflectance  ceases  to  be  independent  of  changes  in 
EQLWC  for  these  cases.  Fig.  9  shows  reflectance  dependence  on  cloud  thickness. 
Channel  1  reflectances  rise  sharply  in  the  first  few  hundred  meters  and  then 
asymptotically  approach  a  value  of  1.0  with  increasing  cloud  thickness.  Channel  3 
reflectances  rise  quickly  at  small  cloud  thicknesses  and  then  remain  constant  beyond  this 
initial  rise.  At  the  point  where  the  reflectance  becomes  constant  with  increasing  thick¬ 
ness,  incident  radiant  energy  has  reached  the  deepest  point  in  the  cloud  layer  where  it 
still  may  be  scattered  back  out  the  cloud  top.  When  channel  3  is  analyzed,  it  must  be 
remembered  that  the  source  of  reflected  radiant  energy  is  confined  to  the  upper-most 
regions  of  the  cloud  layer  (approximately  the  top  50-100  meters).  Also,  reflectance  var¬ 
iations  in  channel  3  will  not  display  the  large  deviations  expected  in  channel  1  which 
result  from  the  complex  changes  in  scattering  optical  depth  over  a  finite  horizontal  dis¬ 
tance  in  a  cloud  field.  For  constant  EQLWC  and  solar  zenith  angle,  reflectance  vari¬ 
ations  in  channel  3  should  be  directly  attributable  to  changes  in  droplet  size  distribution. 
Reflectance  in  both  channels  increases  with  increasing  solar  zenith  angle,  with  depend¬ 
ence  on  the  changing  angle  greater  in  channel  3. 


sr.-rj 
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Fig.  7.  Reflectance  verses  Droplet  Radius  for  z  =  0.25  kjn:  #0  =  45°,  EQLWC  - 
0.4 gjmi  (top)  and  Q.Sg/m3  (bottom).  Dl,  D2,  and  D3  are  the  three  modified 


gamma  distributions. 
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REFLECTANCE  REFLECTANCE 


Fig. 


Reflectance  verses  EQLWC  for  z  =  0.25  km:  d0  =  A5°,  modal  radius  = 
4nm  (top)  and  (bottom).  Di,  D2,  and  D3  are  the  three  modified 


gamma  distributions. 


$ 


m  m  r  •  •  •  ■  r  «  _m  •  •  ■ 

* » o o  > '  v  %•  %  o « ‘  h  • 


27 


REFLECTANCE  REFLECTANCE 


TABLE  5 

DELTA-EDDINGTON  REFLECTANCES  FOR  0.6 3,um 


Distribution  D1 
Distribution  D2 
Distribution  D3 


a  =  2.0  y  =  1.19 
a  =  5.0  y  =  2.4 1 
a  =  5.0  y  =  1.30 


rc  =  S.O  fj.m 


rc  =  8.0m/m 


r,;  =  4.0  ixm  rc  =  4.0,um  rc  =  S.O  jj.m  rc  =  8.0  fj.m 

EQLIVC  =  Q.Sgltn3  EQLIVC  =  0.4g/m3  EQLIVC  =  0.8g/m3  EQLIVC  =  QAgj-n 
Dl  I  D2  I  D3  D 1  I  D2  I  D3  D1  D2  D3  D1  D2  D3 


Reflectances  for  Solar  Zenith  Angle  (0O 


0.912 

0.872 

0.770 

0.558 

0.371 

0.264 

0.091 


0.921 

0.885 

0.794 

0.603 

0.430 

0.269 

0.127 


0.954 

0.933 

0.873 

0.728 

0.561 

0.371 

0.178 


0.959 

0.940 

0.8S6 

0.756 

0.606 

0.430 

0.228 


0.943 

0.917 

0.845 

0.678 

0.500 

0.316 

0.145 


0.S36 

0.770 

0.618 

0.371 

0.211 

0.114 

0.043 


0.912 

O.S73 

0.771 

0.559 

0.371 

0.216 

0.092 


0.892 

0.846 

0.728 

0.500 

0.316 

0.179 

0.071 


0.820 

0.749 

0.589 

0.342 

0.195 

0.100 

0.039 


Reflectances  for  Solar  Zenith  An 


0.949 
0.926 
0.861 
0.711 
0.550 
0.375 
0. 190 


O.S53 

0.794 

0.657 

0.429 

0.269 

0.154 

0.067 


0.921 

0.886 

0.794 

0.604 

0.429 

0.327 

0.127 


0.904 

0.862 

0.756 

0.550 

0.376 

0.229 

0.105 


0.839 

0.775 

0.632 

0.402 

0.243 

0.140 

0.061 


0.900 

0.856 

0.744 

0.522 

0.335 

0.191 

0.077 

gle  (Op) 

0.910 
0.871 
0.771 
0.570 
0.395 
0.243 
0. 1 1 2 


=  30° 

0.879 

0.827 

0.700 

0.464 

0.287 

0.157 

0.064 

=  45° 

0.892 ' 
0.845 
0.731 
0.517 
0.345 
0.206 
0.093 


0.679  0.S17  0.781 
0.589  0.745  0.70(, 
0.399  0.583  0.526 
0.195  0.335  0.283 
0.100  0.191  0.157 
0.050  0.100  0.081 
0.022  0.042  0.035 


0.721  0.836  0.803 
0.632  0.772  0.731 
0.457  0.626  0.574 
0.24S  0.396  0.344 
0.140  0.244  0.205 
0.075  0.137  0.113 
0.031  0.059  0.048 


1 

Rellectanccs  for  Solar  Zenith  Angle  (<h)=  60° 

75  0 

0.933  0.965 

0.957 

0.875  0.933 

0.918 

0.863  0.924  0.908 

0.756 

0.860 

0.833 

500 

0.903  0.949 

0.937 

0.825  0.903 

0  883 

0.809  0.890  0.S69 

0.6S7 

0.806 

0.771 

250 

0.825  0.903 

0.8S2 

0.709  0.825 

0.793 

0.687  0.805  0.772 

0.537 

0.682 

0.638 

loo 

0.663  0.793 

0.755 

0.514  0.664 

0.618 

0.4SS  0.635  0.5S9 

0.34.3 

0.483 

0.435 

.Sit 

0.514  0.665 

0.618 

0.361  0.513 

0.465 

0.343  0.483  0.438 

0.2 1 2 

0.338 

0.296 

25 

0.416  0.513 

0.465 

0.235  0.366 

0.323 

0.212  0.3.38  0.296 

0.122 

0.213 

0.181 

_ 

lo 

0. 198  0.321 

0.278 

O.  IU6  0.199 

0.  i<)2 

0.09X  0.173  0.149 

0.059 

0.103 

0.oS‘ 

/.  /.’  •  ,v«  f.  * 
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TAME  6 

DELTA-EDD1NGT0N  REFLECTANCES  FOR  3.70 pm 


r.  =  4.0  urn 


Distribution  Dl 
Distribution  D2 
Distribution  D3 


r.  =  4.0Mm 


=  2.0  >-=  1.19 
a  =  5.0  y  =  2.41 
a  =  5.0  y  —  1.30 


rc  —  S.O  am 


rc  =  S.Qnm 


EQL (VC  =  O.Sg/m  EOLIVC  =  0.4g/mJ  EQUVC  =  O.Sg/m  EQUVC  =  Q.lg/m 


D3  Dl  D2  D3  Dl  D2  193  Dl 


Reflectances  for  Solar  Zenith  Angle  ((>,,)  =  30° 


D2  D3 


750  0.236 


0.329  0.286 
0.329  0.286 
0.329  0.286 
0.329  0.2S5 
0.316  0.268 
0.258  0.208 
0.144  0.113 


Reflectances  for  Solar  Zenith  Angle  (0„)  =  45' 


.203  0.17 
.203  0.17 
.203  0.17 
.203 


0.169 
0. 1 69 
0.169 
0.169 
0.164 
0.138 
0.083 


0.1 38  0. 10( 


0.176  0.141 
0.128  0.091 
0.069  0.041 


Reflectances  for  Solar  Zenith  .Angle  (10  =  60° 


292  0.428  0.384 
292  0.428  0.384 
292  0.428  0.384 
2S9  0.427  0.383 
270  0.419  0.372 


III.  DATA  ANALYSIS 


A.  SATELLITE  DATA  PROCESSING 

This  section  discusses  the  methods  by  which  the  AVHRR  satellite  data  are  proc¬ 
essed  in  order  to  derive  the  desired  reflective  and  thermal  properties  and  produce  the 
associated  images.  Processing  is  performed  at  the  L'.S.  Naval  Postgraduate  School's 
Interactive  Digital  Environmental  Analysis  (IDEA)  Laboratory.  The  multichannel 
analysis  routine  utilized  in  this  thesis  is  a  FORTRAN  software  package  run  on  a  Digital 
Equipment  Corporation  VAX  8250  computer  system.  The  Gould  DeAnza  image  anal¬ 
ysis  workstation  is  hosted  by  a  Gould  IP8500  image  processor  and  a  9  Mb  micro VAX-I  l 
computer  tied  to  the  VAX  8250  system.  Each  case  study  analysis  in  this  thesis  consists 
of  the  channel  1  radiance  corrected  for  anisotropy,  channel  3  reflectance  corrected  for 
anisotropy,  and  channel  4  brightness  temperature.  Each  analysis  requires  approximately 
45  minutes  processing  tune  (15  minutes  for  channel  1,  25  for  channel  3  and  5  for  channel 
4).  In  comparison,  processing  time  for  any  single  raw  AVHRR  channel  takes  about  5 
minutes.  The  processing  time  includes  reading  the  raw  HRPT  data  from  tape,  calculat¬ 
ing  the  corrected  and  scaled  pixel  values,  and  generating  the  image  for  each  channel.  A 
composite  image  is  constructed  by  overlaying  data  from  the  three  channels.  Each  image 
consists  of  a  512  by  512  pixel,  full  resolution  subscene  located  at  geographic  coordinates 
specified  by  an  initial  satellite  overview  analysis.  Satellite  ephemeris  data  and  subsccne 
quality  are  examined  to  ensure  the  effects  of  sunglint  are  minimized  at  specific  regions 
of  interest. 

1.  Channel  1 

The  measured  radiance  in  channel  1  is  calibrated  in  terms  of  albedo.  The  di¬ 
rectional  reflectance  measured  by  the  satellite  is  converted  to  isotropic  reflectance  by 
weighting  each  pixel  value  by  the  cosine  of  the  solar  zenith  angle  and  the  anisotropic 
reflectance  factor  discussed  in  Chapter  II.  The  resulting  reflectance  values  are  linearly 
scaled  from  20  to  100  percent  for  pixel  values  from  0  to  255.  This  scaling  ensures  that 
there  is  no  rcllcctancc  contribution  from  the  ocean  surface.  The  scaling  also  excludes 
cloudy  pixels  of  minimal  scattering  optical  depth  and  pixels  with  an  insufficient  cloud 
cover  fraction  within  the  field  of  view,  and  enhances  the  contrast  between  stratocumulus 
cloud  tops  of  variable  reflectance. 


Allen  (1987)  utilized  anisotropic  reflectance  factors  in  determining  reflectance  in 
channels  1  and  3  by  selecting  the  appropriate  bin  from  the  data  compiled  by  Taylor  and 
Stowe  (1984)  that  encompassed  the  satellite  zenith  angle,  solar  zenith  angle,  and  relative 
azimuth  for  each  pixel  in  a  subscene.  This  created  problems  at  bin  boundaries  where 
discontinuities  in  anisotropic  reflectance  factors  and  in  the  associated  images  were  ob¬ 
served.  For  channel  1  (and  channel  3)  in  this  thesis,  these  discrete  factors  are  smoothed 
by  entering  them  into  a  three-dimensional  matrix  and  applying  a  moving  least-squares 
fit.  The  resulting  anisotropic  reflectance  factors  are  described  by  basis  functions  whose 
polynomial  coefficients  arc  the  matrix  solution  vector.  The  smoothing  routine  can  in¬ 
corporate  the  factors  designated  for  the  specific  surface  type  being  processed,  in  this  case 
low  level  water  clouds.  Utilizing  the  smoothed  anisotropic  reflectance  factors,  resultant 
isotropic  reflectance  values  are  more  accurate  at  bin  boundaries  by  as  much  as  50  per¬ 
cent  of  the  difference  in  reflectance  between  adjacent  bins  when  no  smoothing  is  applied. 
Discontinuities  in  the  images  due  to  varying  anisotropic  reflectance  factors  arc  elimi¬ 
nated.  There  is  some  concern  that  this  polynomial  method  may  become  unstable  in 
certain  situations;  however,  no  significant  instabilities  arc  observed  in  the  analyses  con¬ 
tained  in  this  thesis.  In  addition,  variations  in  anisotropic  reflectance  as  a  function  of 
wavelength  must  be  ignored  since  the  values  provided  by  Taylor  and  Stowe  (1984)  arc 
spectrally  integrated  values.  The  smoothing  routine  adds  approximately  four  minutes 
to  subscene  processing  time. 

2.  Channel  3 

Radiance  in  channel  3,  comprised  of  both  ;  effected  and  emitted  radiant  energy, 
is  calibrated  in  units  of  mW/m2-sr-arrl.  The  thermally  emitted  portion  is  removed  from 
the  total  radiance  using  a  method  utilized  by  Allen  (1987).  This  method  uses  the 
brightness  temperature  from  channel  4  to  solve  for  the  thermal  emission  in  channel  3. 
Utilizing  Eq.  5,  the  thermal  emission  in  channel  3  is  estimated  by  solving  for  the  Planck 
function  based  upon  the  channel  4  measured  radiance  as  follows: 


/?( T) 


(Cl  14) 


e 


W7/J) 


12) 


The  temperature  is  found  by  manipulating  the  Planck  function  and  is  used  in  Eq.  <i  to 
deduce  the  thermal  emission  in  channel  3.  The  emittance  in  channel  3  is  related  to 


The  fact  that  the  combination  of  thermal  emission  and  reflection  describes  the  total  ra¬ 
diance  in  channel  3  and  that  the  emissivity  in  channel  4  is  close  to  i.O  increases  the  va¬ 
lidity  of  the  estimation  of  the  emissivity  in  channel  3.  Combining  Eq.  5.  10,  and  13 
reveals: 

Lt(cm)  ~  Ach3>£  l  —  P(cmy^  +  WP(C//3)t^(c//3)  cos  ( 14) 

The  reflectance  in  channel  3  is  now  available  utilizing 

_ Ll(CH3)  ~  B(T){C11 3) _ 

“  {,  Ci,cH3)  cos  0„]  -  B(  r)(Cm,  ' 

which  is  the  relationship  used  in  the  analysis  code.  The  anisotropic  reflectance  factors 
in  channel  3  are  the  same  smoothed  function  values  used  in  channel  1  and  discussed  in 
the  previous  section.  Reflectance  values  of  0.0  to  0.45  are  scaled  to  pixel  values  ofO  to 
255  in  the  analysis.  In  deriving  reflectances  from  pixel  values,  a  count  of  0  in  channel 
3  is  assumed  at  any  pixel  where  the  associated  channel  1  reflectance  is  below  20  percent 
(i.e.,  scaled  to  0)  since  the  optically  thin  clouds  in  these  pixels  are  discounted. 

All  channel  3  images  are  slightly  degraded  by  instrument  noise  internal  to  the 
satellite  sensor.  This  noise  appears  in  each  image  as  faint  vertical  bands  of  sinusoidally 
varying  amplitude  in  the  horizontal.  Since  all  of  the  images  analyzed  in  conjunction  with 
this  thesis  are  daytime  passes,  the  effect  of  the  interference  is  minimized. 

3.  Channel  4 

Radiance  values  in  channel  4  are  calibrated  in  units  of  m\V !  m2-sr-cnr' .  Cor¬ 
rections  provided  by  Lauritson  et  al.  (1979)  to  convert  from  radiance  to  temperature  in 
the  NOAA-9/I0  AVIIRR  sensors  are  incorporated  into  the  processing  scheme.  A  tem¬ 
perature  in  the  range  of  200-320  K  is  assigned  to  each  pixel  based  on  the  associated  ra¬ 
diance  value.  The  temperature  values  are  used  to  determine  and  subtract  the 
contribution  of  emitted  radiance  in  channel  3  as  discussed  above.  In  addition,  a 
threshold  temperature  of  273.15  K  is  set  for  image  generation,  so  pixels  that  identify 
below-freezing  temperatures  arc  assigned  a  value  of  255  while  all  others  arc  assigned  a 
value  ofO.  The  resultant  image  highlights  regions  where  cold  clouds  above  the  MAIM, 
may  be  contributing  to  the  radiative  characteristics. 

4.  Composite  Images 

Upon  completion  of  image  generation  for  each  of  the  three  individual  channels, 
a  red-grcen-blue  (RGB)  composite  image  is  produced  by  assigning  channel  3  reflectance 


to  red,  channel  4  threshold  temperature  highlighting  to  green,  and  channel  1  radiance  to 
blue;  and  then  overlaying  the  three  images.  Pixel  values  from  the  individual  analyses  are 
preserved. 

The  resultant  RGB  image  may  provide  a  quick,  general  indication  of  the  nature 
and  variability  of  possible  cloud  droplet  distribution  spectra  within  the  stratocumulus 
clouds.  In  regions  where  the  image  appears  blue,  channel  1  radiance  is  dominating  the 
channel  3  reflectance.  The  relative  reflectance  values  may  be  indicating  a  droplet  si/e 
distribution  composed  of  relatively  large  droplets.  Color  changes  to  violet  and  then  to 
red  indicate  increasing  reflectance  in  channel  3  and,  or  decreasing  reflectance  in  channel 
1.  The  percentage  change  in  reflectance  in  channel  1  for  stratocumulus  cloud  fields  of 
significant  optical  depth  is  much  less  than  the  percentage  change  observed  in  channel 
3;  therefore,  changes  in  channel  1  reflectance  (blue)  due  to  variations  in  optical  depth 
within  the  stratocumulus  cloud  field  are  dominated  by  changes  in  channel  3  reflectance 
(red),  possibly  due  to  variations  in  the  microphysical  properties  of  the  clouds.  Changes 
from  blue  to  red  result  from  a  decrease  in  channel  1  reflectance  due  to  a  loss  of  optical 
depth  and/or  an  increase  in  channel  3  reflectance  which  may  indicate  a  shift  in  the 
droplet  size  distribution  spectra  toward  smaller  droplets.  Green  or  greenish-white  areas 
in  the  RGB  image  indicate  locations  where  cold  clouds  exist.  The  location  of  these  cold 
clouds  are  regions  where  the  scattering  mechanisms  illustrated  in  Chapter  II  become 
more  complex.  In  addition,  these  clouds  are  not  found  in  the  MABL  along  with  the 
stratocumulus  clouds  upon  which  the  analysis  is  focused.  The  highlighted  green  areas 
of  the  composite  image  allow  these  areas  to  be  avoided.  In  addition,  the  ice  cloud 
masking  assists  in  indicating  where  cloud  shadows  arc  modifying  the  reflectance  in  the 
visible  and  NIR.  The  ocean  backgroung  appears  black  with  the  exception  of  isolated 
cases  of  sunglint,  where  the  ocean  appears  as  a  homogenious  red  due  to  the  contribution 
from  channel  3. 

The  major  mechanisms  responsible  for  the  variations  in  droplet  si/c  distribution 
spectra  observed  throughout  the  composite  images  lor  the  case  studies  in  this  thesis  are 
discussed  in  the  next  chapter. 

B.  AIRCRAFT  DATA 

This  section  describes  data  acquired  by  the  NCAR  lilectra  for  missions  coincident 
with  the  case  studies  in  this  thesis.  The  aircraft  data  are  processed  by  NCAR  and  made 
available  in  the  form  of  time-series  plots,  tabuiar  listings,  and  one-minute  interval 
histograms  on  35  mm  microfiche,  accompanied  by  an  applicable  data  quality  report. 


The  data  were  interpreted  using  guidelines  for  standard  output  data  provided  by  Miller 
and  Friesen  (1987).  Aircraft  instruments  and  their  functions  are  as  follows: 

APN-159  Geometric  (Radar)  Altitude. 

Corrected  Particle  Measuring  Systems  (PMS)-King  Probe  EQLWC, 

PMS  Forward  Scattering  Spectrometer  Probe  (FSSP)  Mean  Cloud  Droplet  Diameter, 
PMS  FSSP  Cloud  Droplet  Concentration, 

PMS  Active  Scattering  Aerosol  Spectrometer  Probe  (ASASP)  Mean  Aerosol  Diam¬ 
eter.  and 

PMS  ASASP  Aerosol  Concentration. 

All  PMS  data  in  this  thesis  are  from  one-dimensional  probes  operated  at  a  sampling  rate 
of  one  sample  per  second.  Reflectance  data  used  in  comparison  with  aircraft  measure¬ 
ments  are  from  NOAA-9/10  overpasses  utilized  in  respective  case  studies  and  are  sam¬ 
pled  along  the  approximate  coordinates  of  aircraft  mission  tracks.  The  PMS  data 
displayed  in  Chapter  IV  are  from  representative  postions  along  the  aircraft  sampling 
track,  and  are  interpreted  in  the  context  of  data  in  adjacent  time  intervals  in  order  to 
reduce  the  errors  in  data  interpretation  caused  by  instrument  noise.  It  should  be  noted 
that,  with  the  exception  of  the  14  July  data,  reflectance  cross-sections  coincident  with 
aircraft  data  are  taken  from  a  single  AVHRR  swath,  resulting  in  a  small  difference  in 
azimuth  angle  between  tire  path  of  a  particular  doss-scction  and  the  associated  aircraft 
mission  leg.  Additionally,  there  may  be  several  hours  time  difference  between  aircraft 
sampling  time  and  satellite  overpass.  Due  to  the  priority  in  selecting  case  study  sub¬ 
scenes  (e.g.,  avoiding  broad  areas  of  sunglint)  and  the  difficulty  involved  in  correlating 
reflectance  and  microphysical  data  sampled  at  separate  times  for  individual,  non¬ 
stationary  cloud  elements,  the  18  July  subscene  encompasses  a  region  north  of  the  air¬ 
craft  mission  track;  however,  the  images  show  conditions  analogous  to  those  along  the 
aircraft  mission  track  where  actual  reflectance  sampling  is  performed.  The  PMS  ASASP 
aerosol  probe  was  not  mounted  on  the  aircraft  during  the  29  and  30  June  missions. 


IV.  CASE  STUDY  RESULTS  AND  EVALUATION 


This  chapter  describes  the  four  FIRE  IFO  19S7  case  studies  utilized  in  this  thesis. 
Their  selection  was  made  based  upon  the  availability  of  stratocumulus  clouds  of 
adaquate  horizontal  extent,  sufficient  satellite  coverage,  and  coincidence  with  FIRE 
aircraft  operations.  They  are  dcliniated  as  follows: 

Case  Study  1  (CSI):  2S-30  June 
Case  Study  2  (CS2):  13-14  July 
Case  Study  3  (CS3):  6  July 
Case  Study  4  (CS4):  IS  July 

For  each  case  study,  synoptic/mesoscale  weather  features  are  identified  to  locate  general 
advection  patterns.  Trajectory  analyses  are  included  where  a  more  detailed  How  de¬ 
scription  is  appropriate.  Satellite  images  are  illustrated,  and  cross-sections  of  reflectance 
are  provided  to  assist  in  describing  significant  features.  The  satellite  images  in  this 
chapter  consist  of  channel  1  and  3  reflectances;  refer  to  the  Appendix  for  associated 
RGB  composite  images.  Each  satellite  image  contains  line  segments  or  rectangles  la¬ 
beled  with  single  capital  letters.  These  line  segments  identify  the  location  of  image 
cross-sections  from  which  quantitative  reflectance  data  arc  retrived  and  analyzed.  Air¬ 
craft  data  from  the  NCAR  Electra  are  used  to  relate  cloud  reflectances  obtained  from 
selected  image  cross-sections  to  the  cloud  microphysical  properties  discussed  in 
Chapter  II.  All  times  are  given  in  universal  time  coordinates  (LTC). 

A.  FIRE  IFO  1987  CASE  STUDIES 
1.  CSI  (28-30  June) 

a.  Synoptic  Summary 

The  conditions  off  the  west  coast  of  the  L'.S.  for  the  week  proceeding  CSI 
arc  dominated  by  offshore  flow  over  the  Pacific  Northwest  as  a  result  of  an  inland 
incursion  of  the  subtropical  high.  These  conditions  evolve  toward  a  more  climatological 
northwesterly  flow  during  CSI.  On  28  June,  the  subtropical  high  extends  into  southern 
Alaska  and  is  centered  at  44°N  139° W.  Weak  flow  exists  at  500  mb  over  the  operating 
area.  A  weak  trough  at  500  mb  forms  over  the  northern  California  coast  on  the  29th. 
with  an  associated  weak  (1009  mb)  surface  trough  appearing  as  an  extension  of  the 
persistent  thermal  low,  as  shown  in  Fig.  10.  This  upper-level  trough  intensifies  on  30 
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June  and  moves  to  a  position  centered  at  35°N  128°W,  revealed  in  Fig.  11.  The  center 
of  the  subtropical  high  shifts  westward  to  150°W,  with  an  accompanying  ridge  extending 
towards  Baja  California.  The  surface  trough  merges  with  inland  thermal  troughing  and 
extends  to  Oregon.  The  combination  of  subtropical  high  and  thermal  low  maintained 
westerly  to  northwesterly  surface  winds  in  the  operating  area  over  the  period.  The  San 
Nicolas  Island  rawinsonde  soundings  shown  in  Fig.  12  illustrate  the  maintenance  of  a 
strong  subsidence  inversion  at  the  top  of  the  MABL  and  the  associated  dry  air  at  middle 
and  high  levels. 

Observing  the  synoptic  history  for  this  period,  it  is  obvious  that  there  must 
be  some  offshore  influence  on  the  MABL  due  to  continental  air.  This  is  confirmed  by 
the  NORAPS  trajectory  analyses  displayed  in  Figs.  13  and  14.  The  925  and  S50  mb 
trajectories  near  northern  California  and  Oregon  continue  to  reflect  offshore  flow,  even 
though  the  ridge  driving  the  offshore  flow  is  much  weaker  than  during  the  previous 
week.  The  trajectories  for  30  June,  shown  in  Figs.  15  and  16,  reflect  the  weakening  of 
the  offshore  flow  regime.  They  also  suggest  that  the  continentallv  influenced  air  which 
exists  over  the  ocean  is  being  advected  soutii  and  east  as  the  low-level  winds  become 
more  westerly  over  the  southern  portion  of  the  operating  area. 

The  trajectory  data  in  CS1  (and  CS2)  are  representative  of  the  lowest  layers 
in  the  NORAPS  model.  Nearly  all  of  the  trajectories  reflect  subsidence  as  a  result  of  the 
influence  of  the  subtropical  high  to  the  west  of  the  operating  area.  Since  the  pressure 
changes  for  intermediate  trajectory  locations  include  the  effect  of  the  forecast  model 
levels  changing  with  time,  no  attempt  is  made  to  trace  the  vertical  position  of  the  air 
parcels.  The  pressure  level  given  in  the  applicable  figures  corresponds  to  the  24-hour 
forecast  position.  The  trajectory  data  are  also  compared  to  synoptic  analyses  corre¬ 
sponding  to  the  NORAPS  24-hour  forecast  period  to  ensure  that  NORAPS  forecasts 
verify  with  actual  synoptic  conditions.  The  NO  RAPS  forecasts  are  considered  valid  if 
forecast  winds  are  in  reasonable  agreement  with  verification  analyses. 
b.  AVHRR  Subscenes 

The  satellite  subsccnes  utilized  in  the  CS1  analysis  are  listed  in  Table  7.  The 
coordinates  for  each  subscene  were  selected  during  processing  in  an  attempt  to  highlight 
cloud  microphvsical  gradients,  encompass  NCAR  Llcctra  sampling  legs  where  applica¬ 
ble,  and  avoid  sunglint. 

Figs.  17  and  18  illustrate  channel  1  and  channel  3  reflectance  for  the  28  June 
NOAA-10  satellite  subscene.  The  visible  image  reveals  a  broad,  relatively  homogenious 


Fig.  10. 


29  June  1987  12  UTC  Regional  OI  Analysis:  500  mb  heights  (top)  ami 
surface  pressure/ 1000-500  mb  thickness  (bottom). 


Fig.  11.  30  June  1987  12  UTC  Regional  01  Analysis:  500  mb  heights  (top)  and 

surface  pressure;  1000-500  mb  thickness  (bottom). 


Fig.  13.  NO  RAPS  29  June  1987  00-24  UTC  Parcel  Trajectories:  1000  mb  pres¬ 
sure  at  24  UTC  (top)  and  925  mb  pressure  at  24  UTC  (bottom),  circle 
indicates  24  UTC  position. 


Fig. 


15.  NORAPS  30  June  1987  00-24  UTC  Parcel  Trajectories:  1000  mb  pres¬ 
sure  at  24  L'TC  (top)  and  925  mb  pressure  at  24  LTC  (bottom),  circle 
indicates  24  LTC  position. 
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stratocumulus  cloud  field  over  the  southwestern  half  of  the  subscene.  To  the  northwest 
of  subscene  center,  a  band  of  solid  stratocumulus  is  followed  by  a  band  of  widely  scat¬ 
tered  cumulus.  Finally,  the  northwest  corner  of  the  subscene  is  encompassed  by  another 
solid  stratocumulus  field. 

The  channel  3  image  reveals  several  features  consistent  with  changes  in 
cloud  droplet  size  distribution  within  the  various  cloud  fields.  There  is  a  marked 
boundary  from  southwest  to  northeast  which  appears  to  define  the  extent  of  the  conti¬ 
nental  air  mass  influence.  Based  on  the  reflectance  differences  illustrated  in  Chapter  II, 
the  continental  region  should  appear  significantly  brighter  than  the  marine  regime.  The 
stratocumulus  in  the  immediate  vicinity  of  the  coast,  as  well  as  the  stratus  which  appears 
inland  over  the  Salinas  valley,  is  heavily  influenced  by  continental  CCN  and  has  the 
brightest  reflectance  in  the  image.  In  addition,  a  ship  track  appears  in  the  upper  left 
portion  of  the  image  (coincident  with  cross-section  C).  These  tracks,  as  discussed  in 
detail  by  Coakley  et  al.  (1987),  may  be  the  result  of  the  injection  of  CCN-rich  ship  stack 
exhaust  into  the  cloud  layer.  The  additional  CCN  modify  the  cloud  droplet  size  dis¬ 
tribution,  creating  a  much  higher  concentration  of  smaller  droplets.  As  the  CCN  and 
modified  droplets  are  mixed  to  the  upper  levels  of  the  cloud  layer,  these  regions  can  be 
identified  by  observing  the  higher  reflectance  in  channel  3.  Actually,  the  reflectance  in 
channel  1  also  increases  (as  first  reported  by  Coakley  et  al.  (1987)  and  shown  in  the 
delta-Eddington  model  cloud  reflectances),  but  the  effect  is  usually  masked  by  optical 
depth  variations  in  the  vicinity  of  the  ship  track  resulting  from  EQLWC  and  cloud 
thickness  variations. 


TABLE  7 

CSI  SATELLITE  SL'BSCEN'ES 

DATE 

PLATFORM 

LA  riTL'DE  35°N 
CROSSING  TIME 

SL'BSCENE  CENTER 
COORDINATES 

TIG  ERE 

28  June 

NOAA-10 

16:38:30  ETC 

34"  3 1.LN  1 23°42.6  \V 

mmm 

29  June 

NOAA-9 

22:29:40  ETC 

3 1°  26. ON  1 24°0.7'W 

■ 

30  June 

NOAA-10 

1 5:54:55  ETC 

3 1  "54. 7' N  120°42.4*W 

HBI 

Figs.  19  and  20  show  channel  1  and  channel  3  reflectance,  respectively,  for 
the  29  June  NOAA-9  overpass.  The  location  is  near  the  western  edge  of  the  satellite's 
field  of  view,  causing  the  resolution  of  the  image  on  the  west  edge  to  be  degraded.  The 


NOAA-IO  AVHRR  1639  UTC  28  June  1987  Channel  1:  Line  segments 
indicate  reflectance  cross-sections  A,  B,  and  C.  Center  coordinates 
34°3 1.  l'N  123°42.6'W. 
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Fig.  21.  NOAA-IO  AVHRR  1555  UTC  30  June  1987  Channel  1:  Line  segments 
indicate  reflectance  cross-sections  F,  G,  and  H.  Center  coordinates 
3  l°54.7'N  120°42.4'W. 
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Fig.  22.  NO  A  A- 10  AVHRR  1555  UTC  30  June 

Reflectance:  Line  segments  indicate  reflectance  cross-sections  F.  G.  and 
H.  Center  coordinates  31  °54.7'N  120°42.4'W. 


channel  1  image  reveals  that  the  subscenc  is  mostly  dominated  by  scattered  and  broken, 
optically  thin  stratocumulus  clouds.  The  clouds  appear  slightly  more  broken  in  this 
image  since  some  of  the  optically  thin  clouds  and  partially  cloudy  pixels  of  a  reflectance 
less  than  20  percent  are  scaled  to  zero  by  the  processing  technique.  However,  there  are 
large  clear  areas  in  the  southwest  portion  of  the  subscenc.  The  northwest  corner  of  the 
image  shows  an  encroaching  stratocumulus  field  of  a  more  homogenious  nature.  Point 
Conception  is  visible  in  the  upper  right. 

The  brightness  contrast  in  the  channel  3  image,  in  conjunction  with  the 
continental; marine  interface  observed  in  the  28  June  subscene,  continues  to  be  apparent. 
The  boundary  extends  from  the  left  center  of  the  subscenc  through  the  center  and  up  to 
the  northeast  corner.  The  shift  from  marine  to  continental  regime  is  associated  with  an 
increase  in  the  density  of  cloud  coverage.  The  bright  areas  south  and  west  of  Point 
Conception  are  the  result  of  sunglint  from  the  ocean  surface. 

Images  for  the  30  June  NOAA-10  overpass  appear  in  Figs.  21  and  22. 
Channel  1  reveals  extensive  stratocumulus  cloud  coverage,  becoming  more  broken  to¬ 
wards  the  northwest.  The  cloud  field  in  the  northeast  corner  has  been  modified  by  off¬ 
shore  flow  into  complex  variations  of  optical  depth.  The  offshore  flow  is  a  combination 
of  northwesterly  flow  east  of  Point  Conception,  and  the  morning  extent  of  the  land 
breeze  off  the  Los  Angeles  basin.  The  numerous  dark  spots  over  the  southeast  portion 
of  the  image  are  due  to  cirrus  shadows  (these  can  be  categorized  by  reviewing  the  com¬ 
posite  image  for  30  June  (Fig.  78)  in  the  Appendix). 

The  continental,  marine  boundary  is  still  apparent  in  the  channel  3  image  for 
30  June.  Again,  the  boundary  is  also  associated  with  a  shift  towards  more  dense  cloud 
coverage  in  the  continental  air.  As  expected,  the  region  under  the  influence  of  offshore 
flow,  being  rich  in  continental  CCN,  shows  considerable  reflectance  at  this  wavelength. 
The  bright  spots  in  the  extreme  northeast  corner  of  the  subscene  are  over  the  Los 
Angeles  basin  and  are  most  likely  due  to  industrial  sources  of  CCN  in  convective  plumes 
being  mixed  to  the  top  of  the  cloud  layer.  Another  interesting  feature  in  the  marine  re¬ 
gime  are  the  thin  bands  of  higher  reflectivity  around  individual  cloud  cells.  This  suggests 
support  for  the  interpretations  made  by  Coakley  and  Davies  (1986)  that  increased 
channel  3  reflectivity  at  cloud  edges  should  be  expected  from  the  decrease  in  LQLWC 
and  cloud  droplet  size  distribution  at  the  cloud  edge  where  the  greatest  entrainment  of 
sub-saturated  ambient  air  is  occurring.  This  effect  is  discussed  in  CS4  where  it  is  the 
focus  of  the  case  study. 


Fig.  23  shows  the  progression  of  the  continental/marine  interface  over  the 
3-day  period  constituting  CS1  as  determined  from  satellite  data.  The  geographic  shift 
in  the  boundary  is  consistent  in  direction  with  the  NORAPS  trajectories;  however,  the 
interface  advccts  at  a  slightly  slower  rate  than  the  wind  field  suggesting  that  continental 
aerosol  is  being  freely  mixed  into  the  encroaching  marine  air.  This  also  supports  the 
conclusion  of  Hindman  et  al.  (1977)  that  an  increase  in  the  concentration  of  marine 
CCN  will  not  significantly  affect  the  cloud  droplet  size  distribution  if  a  sufficient  quan¬ 
tity  of  small  (i.e.,  continental)  CCN  remain  in  the  region.  The  additional  possibility  that 
some  continental  aerosols  may  be  entrained  into  the  MABL  from  above  the  subsidence 
inversion  through  turbulent  processes  cannot  be  discounted, 
c.  Reflectance  Cross-Sections 

Channel  1  and  3  reflectance  values  for  cross-sections  A  and  B  are  displayed 
in  Fig.  24.  The  pixel-to-pixel  variability  in  channel  1  is  expected  due  to  the  dependence 
of  reflectivity  in  this  channel  to  variations  in  cloud  thickness.  There  is  a  general  de¬ 
creasing  trend  in  channel  l  reflectance  values  towards  contincntally  influenced  air.  from 
approximately  0.95  to  0.75  over  the  150  pixel  section.  This  effect  is  the  opposite  of  that 
expected  solely  due  to  changes  in  cloud  droplet  size  distribution.  A  trend  towards  de¬ 
creasing  EQLWC  and/or  cloud  thickness  in  the  continental  regime  must  exist  to  account 
for  the  reflectance  variation  in  channel  1. 

The  channel  3  reflectance  values  increase  from  the  marine  to  the  eonti- 
nentally  influenced  air,  from  0.1  to  0.17  for  cross-section  A  and  from  0.1  to  0.23  for 
cross-section  B,  consistent  with  the  anticipated  shift  to  smaller  cloud  droplet  sizes.  Most 
of  the  change  in  reflectance  occurs  over  a  50  pixel  span.  Based  upon  model  reflectances, 
a  change  in  cloud  droplet  modal  radius  of  about  2  to  3^m  exists  across  the 
continental,  marine  interface. 

Cross-section  C,  shown  in  Fig.  25,  highlights  the  reflective  properties  of  a 
ship  track.  The  track  is  not  detectable  in  channel  1  due  to  the  brightness  of  the  sur¬ 
rounding  stratocumulus  cloud  field.  In  channel  3  there  is  a  well-defined  increase  in 
reflectance,  from  an  ambient  value  of  approximately  0.10  to  a  peak  value  of  0.145  at  the 
twelfth  pixel.  Coukley  et  al.  (1987)  report  differences  in  total  channel  3  radiance  (in¬ 
cluding  the  emitted  portion)  of  0.05  to  0.15  m\V,  m2-sr-cm~K  Assuming  no  variability  in 
emission,  these  values  are  consistent  with  those  observed  in  CS1  (as  well  as  the  other 
case  studies  in  this  thesis).  This  difference  represents  a  decrease  in  cloud  droplet  size 
distribution  modal  radius  of  about  1 .5/am  based  on  model  cloud  reflectance  relationships. 


Fig.  23.  CSi  Continental/ Marine  Air  Mass  Interface: 
nei  3  reflectance  from  CSI  subsccnes. 
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Cross-section  D,  displayed  in  Fig.  26,  is  coincident  with  the  region  of  the 
NCAR  Electra  mission  of  29  June.  The  relatively  thin  cloud  field  seen  in  the  satellite 
subscene  is  reflected  in  the  lower  channel  1  reflectance  values  along  the  sampled  path 
(0.33  average).  Channel  3  reflectance  increases  slightly,  from  an  average  of  0.15  to  0.17. 
This  variation  seems  to  indicate  that,  with  decreasing  distance  to  the  landmass,  there  arc 
increases  in  continental  CCN  within  the  air  mass,  but  the  weaker  droplet  size  distrib¬ 
ution  gradient  makes  this  variance  less  noticeable  than  the  obvious  continental,  marine 
air  mass  interface  just  to  the  northwest  of  the  cross-section.  Another  interesting  feature 
is  the  channel  3  reflectance  decrease  associated  with  peaks  in  the  channel  1  reflectance. 
This  may  be  due  to  the  fact  that  the  peaks  in  channel  1  indicate  the  center  of 
stratocumulus  cloud  cells.  Values  of  EQLWC  at  cloud  cell  centers  are  higher  than  the 
values  found  near  cloud  edges  where  subsidence  and  entrainment  mix  sub-saturated  air 
into  the  cloud  from  the  overlying  inversion.  This  mechanism  changes  the  horizontal 
spectrum  of  cloud  droplet  sizes  in  the  upper  50-100  meters  of  the  cloud  (which  is  the 
layer  affecting  the  channel  3  reflectance)  and  creates  the  effect  observed  in  the  figure. 
This  observation  further  supports  the  conclusions  of  Coakley  and  Davies  (I9S6)  con¬ 
cerning  channel  3  reflectance  increases  near  the  cloud  edge  as  a  function  of  decreasing 
cloud  droplet  size  distribution. 

Reflectances  across  a  broken  stratocumulus  cloud  field  are  shown  in  cross- 
section  E,  found  in  Fig.  27.  Again,  there  is  an  obvious  variation  in  channel  3  reflectance 
towards  the  edges  of  the  cloud  fields.  As  the  cross-section  approaches  a  cloud  edge,  the 
channel  3  reflectance  increases  in  conjunction  with  decreasing  droplet  size.  There  is 
usually  a  slight  decrease  in  reflectivity  at  the  immediate  edge,  since  the  clouds  at  this 
point  are  reaching  a  sufficiently  thin  optical  depth  and  no  longer  appear  infinitely  thick 
to  channel  3.  Further  decreases  in  cloud  thickness  cause  decreases  in  channel  3  due  to 
the  diminishing  number  of  cloud  droplets  available  to  scatter  the  incident  radiation. 
Recall  that  the  image  processing  technique  drives  tiie  channel  3  reflectance  value  to  zero 
where  the  associated  channel  1  value  has  been  scaled  to  zero,  in  reality,  there  may  be 
a  small  contribution  to  channel  3  between  the  cloud  regions,  but  it  is  from  the  ocean 
surface  and  hinders  the  ability  to  analyze  the  clouds. 

Satellite-aircraft  geographic  coincidence  occurs  again  on  30  June.  The 
rellectancc  is  shown  in  Fig.  2S  as  cross-section  i:.  Channel  1  decreases  towards  the 
continental  air  to  the  cast,  from  an  average  reflectance  of  0.85  to  0.80.  Channel  3 


wuuummui  m  a  w  wlsw.vwv  v  r?  lt7 


■V  •  T7T7n^T^f7yyy^y»3TO?y  TvTTTO'A'.m  ••wrv.'  <- 


Fig.  26.  Reflectance  Cross-Section  D:  Sampled  west  to  east;  coincides  with  lo¬ 
cation  of  29  June  NCAR  Electra  mission. 


reflectance  increases  from  0.08  to  0.14  across  the  sampled  path,  representing  a  change 
in  droplet  size  distribution  modal  radius  of  approximately  Lo^m. 

Cross-section  G  is  illustrated  in  Fig.  29  and  is  a  good  example  of  the  dif¬ 
ferent  effect  in  the  two  channels  resulting  from  cloud  optical  depth  variations.  Large 
excursions  in  channel  1  reflectance  result  from  the  sensitivity  of  this  channel  to  vari¬ 
ations  in  cloud  thickness.  The  associated  channel  3  values  remain  relatively  immune  to 
these  variations,  but  appear  slightly  more  variable  than  previous  examples  in  conjunc¬ 
tion  with  cloud  microphysical  fluctuations  in  the  horizontal.  Channel  3  reflectance  in¬ 
creases  from  an  average  of  0.08  to  0.16  from  north  to  south  along  the  cross-section. 

The  final  reflectance  cross-section  for  CS1  is  labeled  H  and  is  shown  in  Fig. 
30.  This  track  illustrates  the  effects  of  offshore  flow  in  the  immediate  vicinity  of  the  Los 
Angeles  basin.  The  influence  of  a  substantial  source  ofCCN  of  urban,  industrial  origin 
is  readily  apparent  in  both  channels.  The  channel  1  profile  indicates  a  0.2  drop  in 
reflectance  at  the  extent  of  the  strong  continental  influence  from  the  loss  of  optical  depth 
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associated  with  air  of  less  EQLWC.  Channel  3  rises  sharply  at  the  interface,  from  an 
average  reflectance  of  0.08  to  0.2.  More  complex  circulation  patterns  within  this  region 
of  influence  and  the  accompanying  variations  in  cloud  microphysical  properties  add 
considerable  structure  to  the  channel  3  reflectance.  The  reflectance  gradient  suggests  a 
change  in  cloud  droplet  modal  radius  of  approximately  3 /mi. 
d.  Aircraft  Data 

( I)  29  June  1987.  The  mission  track  for  the  NCAR  Electra  on  29  June 
is  shown  in  Fig.  31,  with  the  track  segment  from  point  C  to  point  B  corresponding  with 
satellite  reflectance  cross-section  D  from  Fig.  26.  The  Electra  data  are  taken  from 
2115-2120  UTC  while  the  aircraft  maintained  a  posture  just  below  the  cloud  top,  as¬ 
cending  from  1040  meters  at  time  2115  L'TC  to  the  approximate  location  of  the  cloud 
top  (1065  meters)  at  time  2120  UTC.  This  corresponds  to  the  center  half  of  the  point 
C  to  point  B  mission  leg.  PMS  data  for  29  June  are  displayed  in  Table  8.  Although  the 
aircraft  track  does  not  cross  any  significant  continental  marine  gradient,  a  general  trend 
towards  higher  droplet  concentration  and  smaller  cloud  droplet  mean  radius  is  observed 
as  the  aircraft  passes  from  west  to  east.  This  is  consistent  with  the  gradual  increase  in 
reflectance  illustrated  in  cross-section  D. 


TABLE  8 

SELECTED  NCAR  ELECTRA  DATA  FOR  29  JUNE  1987 


Time 

(UTC) 

(±20sec) 

APN-159 

Altitude 

(m) 

PMS- King 
EQLWC 
(glnO) 

AS  ASP 
cone. 

(XI  cm3) 

ASASP 

mean 

radius 

(Mm) 

FSSP 

cone. 

(Xian3) 

FSSP 

mean 

radius 

(am) 

■m 

N  A 

N  A 

48 

9.3 

2116 

1035 

N.'A 

N.'A 

54 

8.S 

2117 

1040 

0,30 

N.'A 

N  A 

55 

8.5 

2118 

1055 

0.29 

N  A 

NT  A 

68 

S.4 

2119 

1050 

0.29 

N  A 

N  A 

58 

8.5 

f2j  30  June  1987.  The  NCAR  Electra  mission  track  on  30  June  is  illus¬ 
trated  in  Fig.  32;  the  track  segment  from  B  to  A  corresponding  with  reflectance  cross- 
section  F  found  in  Fig.  28.  The  aircraft  data  are  taken  from  two  cloud  sounding  passes, 
from  2032-2042  UTC  (B  to  A)  and  from  2306-2316  UTC  (A  to  B).  In  these  sounding 
runs,  the  aircraft  "porpoises"  between  cloud  base  and  the  clear  air  in  the  inversion  above 


Fig.  31.  NCAR  Electra  Flight  Track,  29  June  1987:  (from  Kloesel  ei  at.,  1988). 


the  cloud  top.  This  flight  profile  is  especially  useful  in  defining  the  cloud  microphysical 
structure  near  the  cloud  top,  being  the  source  of  channel  3  reflectance;  however,  it  limits 
the  data  describing  horizontal  variations  at  any  one  level.  In  the  first  sounding  leg,  the 
aircraft  climbs  through  the  cloud  top  at  20:33:20,  20:38:30  and  at  20:41:05  UTC;  de¬ 
scends  into  the  cloud  top  at  20:35:40  and  at  20:38:30  UTC,  and  dips  slightly  below  the 
cloud  base  at  20:36:30  and  at  20:40:20  UTC.  In  the  second  sounding  run,  the  aircral't 
exits  cloud  top  near  23:07:50  and  23:13:00  UTC;  descends  into  the  cloud  at  23:10:15  and 
23:15:15  LTC;  and  is  at  the  cloud  base  at  23:06:55,  23:11:55,  and  23:17:05  UTC.  PMS 
data  for  both  sounding  legs  are  shown  in  Table  9. 

In  both  aircraft  sampling  legs,  there  are  indications  of  a  shift  toward 
higher  droplet  concentrations  and  smaller  droplet  radius  at  the  eastern  end  of  the  track. 


This  is  consistent  with  the  trend  expected  from  a  review  of  cross-section  F.  Fig.  33 
shows  droplet  size  distribution  histograms  at  one  sample  per  second  averaged  over  one 
minute  for  the  cloud  samples  closest  to  the  aircraft  track  endpoints  in  both  aircraft 
sampling  runs.  In  both  cases,  the  distribution  nearer  to  point  A  is  shifted  towards  a 
smaller  droplet  modal  radius  and  the  concentration  of  large  droplets  decreases. 

2.  CS2  (13-14  July) 

a.  Synoptic  Summary 

A  general  transition  from  trough  to  ridge  at  the  level  of  non-divergence 
characterizes  the  synoptic  evolution  during  CS2.  On  12  July  a  500  mb  trough  exists  over 
the  Rocky  Mountains  and  continues  its  movement  to  the  east  through  the  period.  A 
ridge  begins  to  gradually  build  at  136°W,  impinging  on  the  coast  near  Washington  and 
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FSSP  Probe  Droplet  Size  Distribution  Histograms  for  30  June 
1987:  First  leg  (top)  and  second  leg  (bottom).  Dashed  lines  are  cloud 
samples  near  point  B,  solid  lines  represent  point  A. 


i  »  a  V  •\>rn 


TABLE  9 

SELECTED  NCAR  ELECTRA  DATA  TOR  30  JUNE  1987 


Time 

(ETC) 

(±20sec) 


APN-159 

Altitude 

(m) 


P.V1  S- King 
EQLVVC 
(. gl mi) 


ASASP 

cone. 

( XI  cm3) 

ASASP 

mean 

radius 

(urn) 

FSSP 

cone. 

(A’/cm') 

N/A 

N/A 

30 

N/A 

N/A 

30 

N/A 

N/A 

38 

N/A 

N/A 

55 

N/A 

N/A 

33 

N/A 

N/A 

37 

N/A 

N/A 

39 

N/A 

N/A 

32 

N/A 

N/A 

30 

FSSP 
mean 
radius 
(t-  tm) 


0.5 

0.0 

1.0 

0.0 

9.0 


British  Columbia.  The  center  of  the  subtropical  high  is  located  at  39°N  142°W,  with  a 
lobe  of  the  high  reaching  the  IFO  operating  area.  The  1010  mb  central  pressure  of  the 
thermal  trough  is  several  millibars  higher  than  the  previous  day,  resulting  in  light, 
northwesterly  surface  winds  over  the  region  of  interest.  Significant  offshore  (low  exists 
from  northern  California  to  Washington  that  decays  in  strength  as  the  day  progresses. 
The  500  mb  ridge  builds  to  5940  meters  on  13  July  and  moves  over  the  coast,  as  the 
surface  high  builds  in  toward  Canada  as  shown  in  Fig.  34.  Also,  the  thermal  trough 
strengthens  to  1009  mb  causing  an  increase  in  the  strength  of  northwesterly  winds  olT 
the  southern  California  coast.  Continued  strengthening  of  the  subtropical  high  and  the 
thermal  trough  (to  1007  mb)  on  14  July,  illustrated  in  Fig.  35,  results  in  strong 
northwesterly  winds  over  the  operating  area.  The  500  mb  pressure  gradients  are  quite 
weak.  Ruwinsondc  soundings  taken  from  the  research  vessel  POINT  SL  R,  operating 
just  upwind  from  San  Nicolas  Island,  are  displayed  in  Fig.  36.  The  soundings  confirm 
the  presence  of  strong  subsidence  and  vertical  decoupling  across  the  stable  inversion  al¬ 
lowing  the  shallow  thermal  trough  to  build  the  surface  winds. 

In  a  situation  similar  to  CS1,  the  NORAPS  trajectory  analyses  shown  in 
Figs.  37  through  40  reveal  a  source  of  continental  air  being  adveeted  into  the  IFO  op¬ 
erating  area.  Strong  offshore  How  exists  over  northern  California  on  12  July  at  925  mb 
up  through  and  including  the  700  mb  level.  A  shear  zone  exists  near  125°W  at  all  levels, 
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but  especially  at  850  and  700  mb  where  a  weak  trough  is  driving  southerly  How  near  the 
coast.  As  the  surface  pressure  gradient  builds  on  13  July,  strong  northwesterly  How  be¬ 
comes  apparent  at  1000  mo,  advecting  air  of  a  continental  nature  into  the  operating 
area.  This  shift  in  the  flow  pattern  is  consistent  with  the  replacement  of  continentaliy 
influenced  air  off  the  northern  California  coast  with  marine  air  from  the  northwest. 
Offshore  winds  are  still  providing  a  source  of  continental  CCN  both  above  and  belowr 
the  subsidence  inversion  over  northern  California  in  the  early  stages  of  the  case  study 
period.  The  evolving  flow  pattern  is  of  a  nature  that  should  confine  continentaliy  in¬ 
fluenced  air  nearer  to  the  coastline  with  decreasing  latitude. 
b.  A  VHRR  Subscenes 

Satellite  s  lbscenes  found  in  CS2  are  listed  in  Table  10.  As  in  CS1,  the  pri¬ 
ority  for  subscene  location  is  the  location  of  cloud  field  reflectance  gradients  indicating 
differences  in  air  mass  qualities. 


TABLE  10 

CS2 SATELLITE  SUBSCENES 

DATE 

PLATFORM 

LATITUDE  35°N 
CROSSING  TIME 

SUBSCENE  CENTER 
COORDINATES 

FIGURE 

13  Julv 

NOAA-10 

16:12:55  UTC 

32°48.7'N  120°29.7'W 

41/42 

13  Julv 

NOAA-9 

23:20:35  UTC 

33°57.S'N  123°30.1'W 

43  44 

14  Julv 

NOAA-10 

15:51:00  UTC 

34°59.8'N  1 24°  16.9'W 

45  46 

14  Julv 

NOAA-9 

23:09:40  UTC 

32°33.7'N  121°48.8'W 

47  48 

The  channel  1  and  3  reflectance  patterns  for  the  morning  of  13  July  are  il¬ 
lustrated  in  Figs.  41  and  42.  The  visible  image  indicates  mixing  of  the  dry  air  above  the 
inversion  into  the  MABL  downwind  of  the  coastal  boundary  cast  of  Point  Conception 
and  the  Channel  Islands  due  to  the  strong  northwesterly  surface  winds.  The 
stratocumulus  cells  are  aligned  in  streets  parallel  to  the  llow,  with  larger  individual  colls 
apparent  in  the  western  extreme  of  the  image.  The  stratocumuius  field  is  slightly  less 
rellective  near  the  coast  north  of  Point  Conception.  The  channel  3  image  shows  an 
unmistakable  contrast  between  the  continental  and  marine  air  masses.  The  brighter  re¬ 
gion  is  broader  to  the  north  where  it  is  closer  to  the  origin  of  the  offshore  llow,  and 
confined  near  the  coast  by  the  advection  pattern  further  to  the  south.  Fvcn  brighter 
regions  are  found  against  the  coast  where  the  continental  CCN  influence  is  greatest. 


Fig.  38.  NORAPS  12  July  1987  00-24  UTC  Parcel  Trajectories:  850  mb  pressure 
at  24  UTC  (top)  and  700  mb  pressure  at  24  UTC  (bottom),  circle  indicates 
24  UTC  position. 


Fig.  39. 


NORAPS  13  July  1987  00-24  UTC  Parcel  Trajectories:  1000  mb  pres¬ 
sure  at  24  UTC  (top)  and  925  mb  pressure  at  24  U'TC  (bottom),  circle 
indicates  24  UTC  position. 


NOAA-IO  AVHRR  1613  UTC  13  July  1987  Channel  1:  Line  segments 
indicate  reflectance  cross-sections  M  and  N.  Center  coordinates 
32°48.7'N  120°29.7'W. 


Fig.  42. 


NOAA-IO  AVHRR  1613  UTC  July  1987  Channel  3  Reflectance:  Line 
segments  indicate  reflectance  cross-sections  M  and  N.  Center  coordinates 
32°48.7'N  120°29.7'W. 


The  subscene  images  from  the  NOAA-9  afternoon  pass,  shown  in  Figs.  43 
and  44,  are  centered  further  north  in  the  operating  area  off  the  south-central  California 
coast.  Channel  1  reveals  a  relatively  uniform  stratocumulus  cloud  field  over  the  region. 
The  clouds  in  the  northwest  half  of  the  image  show  more  of  a  cellular  structure,  with  the 
cell  centers  appearing  the  brightest.  The  clouds  closer  to  land  are  also  structured,  but 
are  more  uniform  in  brightness.  The  high  reflectance  in  channel  3  near  the  coast  is  a 
result  of  the  subscene  being  near  the  source  of  offshore  flow  on  13  July.  The  center  of 
the  image  displays  a  north-south  rellectance  band  which  is  not  as  bright  as  the  air  nearer 
to  the  coast  under  the  immediate  influence  of  continental  aerosols,  but  is  brighter  than 
the  marine  air  which  appears  further  to  the  west.  The  clouds  of  intermediate  brightness 
in  this  image  are  most  likely  the  result  of  air  which  was  advected  offshore  further  north. 
The  clouds  within  this  air  mass  may  be  under  the  influence  of  CCN  of  different  compo¬ 
sition  and  concentration  than  the  air  mass  closer  to  the  coast.  In  addition,  microphys¬ 
ical  modifications  to  CCN  content  during  the  time  the  air  mass  is  advected  over  the 
ocean  may  create  a  small  difference  in  radiative  properties  by  the  time  these  parcels 
reach  the  operating  area.  As  discussed  by  Noonkester  (1984),  continental  air  existing 
over  the  ocean  is  continually  being  modified  by  the  aerosol  loss  processes  described  in 
Chapter  II.  Note  that  this  region  of  intermediate  brightness  has  a  similar  channel  3 
reflectance  to  the  continental  regime  obser'-fd  further  south  during  the  morning. 

The  channel  1  it  age  for  the  morning  of  14  July,  found  in  Fig.  45,  shows  a 
continuous  and  uniform  stratocumulus  cloud  field,  slightly  less  bright  in  the  eastern  third 
of  the  subset  e.  The  eddies  which  appear  in  the  cloud  topography  over  Monterey  Bay 
and  to  the  south  along  the  coast  indicate  that  surface  winds  are  still  relatively  strong, 
and  are  in  a  direction  along  the  coast.  There  are  several  interesting  features  in  the 
channel  3  image,  displayed  in  Fig.  46.  The  most  prominent  is  the  lower  brightness  near 
the  coast  within  the  region  of  continental  influence.  The  clouds  near  the  coast  are 
brighter  than  the  marine  air  to  the  southwest,  as  expected;  but  are  not  as  bright  as  the 
clouds  near  the  center  of  the  image.  The  strong  northwesterly  flow  has  replaced  the 
continental  air  mass  near  the  coast  with  a  surge  of  marine  air  from  the  northwest.  The 
sequence  of  images  seems  to  indicate  that  the  region  of  high  channel  3  reflectivity  ap¬ 
parent  in  the  13  July  NOAA-9  overpass  has  been  advected  rapidly  southward  by  the 
strong  ambient  flow.  The  brighter  region  near  the  center  of  the  image  is  most  likely  the 
result  of  a  supply  of  continental  air  advected  offshore  over  Oregon,  Washington,  and 
Canada. 


Fig.  43. 


NOAA-9  AVHRR  2321  UTC  13  July  1987  Channel  1:  Line  segments 
indicate  reflectance  cross-sections  O  and  P.  Center  coordinates 
33#57.8'N  123°30.rW. 
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Fig.  46.  NOAA-IO  AVHRR  1551  UTC  14  July  1987  Channel  3  Reflectance:  Line 
segments  indicate  reflectance  cross-sections  Q  and  R.  Center  coordinates 
34°59.8'N  124°16.9'W. 
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The  continental  marine  interface  observed  on  13  July  is  along  a  north-south 
line  left  of  the  subscene  center.  The  smaller  cloud  droplet  size  distribution  associated 
with  the  thin  stratus  of  the  Salinas  valley  and  within  the  small  eddy  south  of  Monterey 
are  characteristically  brighter  than  surrounding  clouds.  In  addition,  two  ship  tracks  are 
highlighted,  one  in  the  northwest  corner  of  the  subsccne  and  one  near  the  subscenc 
center. 

Figs.  47  and  48  display  channel  1  and  3  reflectance  for  the  NOAA-9  after¬ 
noon  pass  on  14  July.  The  visible  channel  reveals  the  maintenance  of  continuous 
stratocuinulus  clouds,  with  the  exception  of  the  region  of  the  Channel  Islands.  'Ihe 
southwestern  third  of  the  image  appears  more  cellular  in  nature  than  the  cloud  field 
closer  to  the  continent  with  the  stratocumuius  becoming  more  broken  in  the  extreme 
southwest.  A  mesoscale  gravity  wave  appears  as  a  line  of  discontinuity  along  the  coast 
north  of  Point  Conception.  This  wave  can  be  tracked  using  geostationary  satellite  data 
and  is  observed  to  propagate  northwest  to  its  approximate  position  in  Fig.  47/48  and 
then  dissipates  in  the  ambient  flow.  Low  channel  1  reflectance  due  to  lower  EQLWC 
and  optical  depth  caused  by  entrainment  is  readily  apparent  to  the  southeast  of  this 
feature.  The  channel  3  image  shows  several  regions  of  variable  reflectance.  The  cloud 
field  is  most  bright  south  of  Point  Conception  and  decreases  in  brightness  towards  the 
northwest  across  three  diffuse  boundaries  of  reflectance  magnitude.  A  bright  band  ex¬ 
tends  northwest  to  southeast  in  the  lower  left  quadrant  of  the  subscene,  with  dark  (ma¬ 
rine  influenced)  reflectance  to  the  south  of  this  band.  This  feature  is  intersected  by  the 
NCAR  Electra  mission  track  for  14  July.  As  expected,  the  region  southeast  of  the  wave 
feature  is  brighter  than  the  air  mass  north  of  the  line  of  discontinuity. 

The  continental; marine  interfaces  in  CS2  fluctuate  longitudinally  in  relation 
to  the  wind  field,  as  shown  in  Fig.  49.  The  history  of  continental  surges  of  air  is  more 
complex  in  this  case  study  since  the  surface  winds  are  relatively  strong  and  the  surges 
of  continental  air  are  advected  rapidly  to  the  south.  The  positions  of  the 
continental/marine  interfaces  in  Fig,  49  suggest  that  the  strength  of  offshore  flow  and  its 
north-south  position  along  the  coast  determine  the  westward  extent  and  the  persistence 
of  the  continental  influence. 

c.  Reflectance  Cross-Sections 

Channel  1  and  3  reflectance  values  for  cross-sections  \1  and  N  are  displayed 
in  Fig.  50.  Again,  the  characteristic  trend  towards  tower  reflectance  in  continental  air 
is  seen  in  channel  1,  with  values  decreasing  from  approximately  0.75-0.8  in  the  marine 


7~zrvrrr'<rKr*. 


NOAA-9  AVHRR  2310  UTC  14  July  1987  Channel  1:  Channel  1  (top) 
and  channel  3  reflectance  (bottom),  with  reflectance  cross-sections  S,  T, 
and  U.  Center  coordinates  32°33.7’N  121°48.8'W. 
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NOAA-9  AVHRR  2310  UTC  14  July  1987  Channel  3  Reflectance:  Line 
segments  indicate  reflectance  cross-sections  S,  T,  and  L\  Center  coordi¬ 
nates  32°33.7'N  121°48.8'W. 


Fig.  48. 


Fig.  49. 


CS2  Continental/ Marine  Air  Mass  Interface:  Identified  utilizing  channel 
3  reflectance  from  CS2  subscenes.  Dashed  lines  indicate  diffuse 
reflectance  boundaries. 
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air  to  0.7  across  the  boundary.  The  trend  in  channel  3  in  the  same  direction  is  an  in¬ 
crease  from  around  0.15  to  0.2,  constituting  a  change  in  cloud  droplet  modal  radius  of 
approximately  1.5  to  2^m  based  on  modal  reflectance  calculations.  The  channel  3 
reflectance  is  somewhat  more  variable  from  pixel  to  pixel  in  this  case,  most  likely  due  to 
droplet  size  distribution  variations  as  a  result  of  turbulent  mixing  from  the  strong  surface 
wind  conditions. 

Fig.  51  shows  cross-sections  O  and  P  from  the  13  July  NOAA-9  overpass. 
The  decrease  in  channel  l  in  the  continental  regime  is  again  observed,  and  is  much  more 
evident  in  this  figure.  Channel  1  reflectance  decreases  from  around  0.8  to  0.6  for  cross- 
section  O  and  from  0.75  to  0.5  for  cross-section  P.  The  gradient  in  channel  3  reflectance 
is  also  quite  strong  in  this  case,  increasing  approximately  0.2  for  each  case.  The  inter¬ 
mediate  area  of  reflectance  discussed  in  the  previous  section  has  a  value  near  0.21,  sim¬ 
ilar  to  the  reflectances  obtained  in  the  continental  air  mass  from  the  morning  NOAA-10 
overpass.  The  area  of  "fresh"  continental  CCN  near  the  coast  displays  some  of  the 
highest  channel  3  reflectances  (and  the  strongest  trend  toward  low  channel  1  reflectance) 
in  this  thesis.  It  can  be  expected  that  the  cloud  droplet  size  distribution  within  this  re¬ 
gion  is  strongly  shifted  toward  the  small  end  of  the  anticipated  size  spectrum. 

Cross-sections  Q  and  R,  found  in  Fig.  52,  show  the  conditions  on  the 
morning  of  14  July.  The  channel  1  profiles  show  no  significant  trend,  although  they 
appear  somewhat  less  variable  due  to  the  uniform  nature  of  the  stratocumulus  cloud 
field.  Cross-section  Q  in  channel  3  shows  a  rise  in  reflectance  from  around  0.19  in  the 
marine  environment  to  0.25  in  continentaliy  influenced  conditions.  The  ship  track  at 
pixel  13  is  clearly  visible.  Channel  3  cross-section  R  shows  the  decrease  in  reflectivity 
towards  the  coast,  the  source  of  which  being  the  advection  of  marine  air  into  the  near¬ 
coastal  environment  due  to  the  strong  surface  winds. 

Cross-section  S,  shown  in  Fig.  53,  displays  the  change  in  reflectance  across 
a  diffuse  boundary  of  channel  3  reflectance  intensity.  Channel  1  reflectance  is  typically 
variable  and  remains  near  an  average  value  of  0.7  over  the  cross-section.  The  associated 
channel  3  reflectance  change  is  approximately  0.1  to  0.18  from  marine  to  continental 
influence,  corresponding  to  a  model  cloud  droplet  size  distribution  modal  radius  change 
of  about  2^m. 

Fig.  54  displays  a  short,  25  pixel  cross-section  labeled  as  T,  which  traverses 
the  wave  feature  seen  in  the  14  July  NOAA-9  satellite  data.  The  sharp  decrease  in 
channel  1  in  conjunction  with  decreases  in  EQLWC  and  cloud  optical  depth  dominate 
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Fig.  52.  Reflectance  Cross-Sections  Q  and  R:  Sampled  west  to  east;  ST  denotes 
position  of  ship  track. 
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Fig.  54.  Reflectance  Cross-Section  T:  Sampled  north  to  south. 

this  profile.  Channel  3  reflectance  rises  due  to  the  entrainment  of  sub-saturated  air,  from 
0.12  to  0.2  across  the  boundary  towards  the  south. 

Reflectance  cross-section  U  coincides  with  the  location  of  the  14  July 
NCAR  Electra  mission  and  is  shown  in  Fig.  55.  This  cross-section  intersects  a  moderate 
channel  3  reflectance  gradient  (refer  to  Fig.  48).  The  cellular  nature  of  the  stratocumulus 
clouds  in  this  region  is  graphically  apparent  in  the  channel  I  profile  which  fluctuates  up 
to  30  percent  within  the  sampled  path.  Channel  3  increases  from  an  average  value  of 
0.08  in  the  marine  environment  to  approximately  0.12  across  the  reflectance  boundary, 
representing  a  change  in  cloud  droplet  modal  radius  of  approximately  1-1.5/im  based  on 
model  cloud  reflectance  data.  The  gradient  is  within  the  first  30  pixels  of  the  cross- 
section,  with  a  slight  decrease  around  the  40th  pixel. 
d.  Aircraft  Data 

(])  14  July  1987.  Fig.  56  shows  the  NCAR  Electra  flight  track  for  the 

14  July  cloud  sampling  mission.  This  track  corresponds  to  reflectance  cross-section  U, 
displayed  in  Fig.  55.  The  aircraft  data  are  taken  in  a  cloud  sampling  leg  from  2127-2135 


Fig.  55.  Reflectance  Cross-Section  U:  Sampled  west  to  east;  coincides  with  lo¬ 
cation  of  14  July  NCAR  Electra  mission. 


L'TC  along  a  path  from  point  K  to  point  J.  The  mission  was  run  in  a  solid 
stratocumulus  cloud  layer  with  the  cloud  base  at  640  meters  and  the  cloud  top  at  an 
average  height  of  830  meters.  The  Electra  remained  in  the  upper  half  of  the  cloud  layer 
during  the  sampling  run,  exiting  and  reentering  the  cloud  top  eleven  times  (one  minute 
intervals  starting  at  21:27:20  ETC  for  three  minutes,  then  eight  excursions  at  approxi¬ 
mately  35  second  intervals)  during  the  eight  minute  leg.  Table  11  contains  the  aircraft 
microphysical  data  for  this  leg. 

The  aircraft  measurements  reveal  a  striking  gradient  in  aerosol  con¬ 
centration,  from  1550  to  155  per  cm\  The  mean  aerosol  radius  increases  in  size  to  the 
west  by  0.0 l^m.  Due  to  the  higher  concentration  of  CCN  at  point  K,  the  cloud  droplet 
concentration  is  two  and  one-halftimes  as  large  (85  per  cm?)  as  the  concentration  at 
point  J  (34  per  cm').  The  associated  mean  droplet  radius  increases  toward  the  marine 
environment  by  1.7/atn  within  the  track.  This  compares  favorably  with  the  reflectance 
gradient  in  cross-section  U.  Fig.  57  shows  a  comparison  of  the  aircraft  track  endpoint 
histograms  for  cloud  droplet  and  aerosol  size  distribution.  T  he  histograms  reveal  the 
change  in  concentration  for  both  droplets  and  aerosols  across  the  continental,  marine 
interface.  The  FSSP  histogram  representing  the  continental  environment  shows  a  higher 
concentration  of  small  droplets  and  a  smaller  modal  radius  than  in  the  marine  regime. 
The  ASASP  histograms  show  the  significant  increase  in  small  aerosol  particles  in  the 
continental  air  at  point  K. 


TABLE  1! 

SELECTED  NCAR  ELECTRA  DATA  FOR  14  JULY  1987 


Time 

(ETC) 

(±20sec) 

APN-159 

Altitude 

(m) 

PMS-King 
EQLWC 
(gin P) 

ASASP 

cone. 

(A/ cm3) 

ASASP 

mean 

radius 

(/mi) 

FSSP 

cone. 

(Stem') 

FSSP 

mean 

radius 

(um) 

2127 

805 

msm I 

1550 

0.105 

9.9 

2128 

815 

1050 

IS11BI 

10.8 

2130 

795 

0.30 

450 

jamli  ■ 

44 

9.7 

2131 

825 

0.34 

270 

61 

9.8 

2133 

825 

0.34 

200 

0.101 

58 

10.5 

2134 

820 

0.30 

175 

40 

11.0 

2135 

825 

0.32 

155 

B3E3sB 

34 

11.4 

3.  CS3  (6  July) 

a.  Synoptic  Summary 

The  500  mb  pattern  on  6  July,  as  shown  in  Fig.  58,  is  characterized  by  a 
weak  wave  over  California  imbedded  in  basically  zonal  flow.  The  subtropical  high 
dominates  the  eastern  Pacific  with  ridges  extending  into  the  northwestern  L'.S.  and  into 
the  IFO  operating  area.  The  thermal  trough  is  at  1007  mb,  producing  moderate  winds 
out  of  the  north  over  the  CS3  region.  A  phenomena  known  as  the  Catalina  Eddy  (so 
called  because  it  usually  appears  centered  near  Santa  Catalina  Island)  persists  from  its 
formation  on  the  previous  day.  As  discussed  by  Wakimoto  (1987),  the  cyclonic  shear 
vorticity  which  appears  downstream  of  the  coastal  mountains  leads  to  subsidence  south 
of  Point  Conception  and  Santa  Barbara,  wliile  the  southerly  flow  west  of  San  Diego 
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experiences  convergence  and  lifting.  The  eddy  has  the  tendency  to  ventilate  the  air  in 
the  urban  environment  over  the  Los  Angeles  basin.  Fig.  59  shows  the  1646  ITC  San 
Nicolas  Island  rawinsonde  sounding  and  confirms  the  conditions  expected  under  the  in¬ 
fluence  of  subtropical  high  subsidence  and  the  thermal  trough. 
b.  A  VHRR  Subsce/ie 

A  single  satellite  subscene  is  generated  for  CS3  on  6  July;  the  images  arc 
displayed  in  Figs.  60  and  61.  The  NOAA-lO  satellite  subscene  has  center  coordinates 
of  32°35.2'N  122°33.4  W  and  crosses  latitude  35°N  at  15:24:20  ETC  on  a  decending 
node.  The  channel  1  image  clearly  shows  a  spiral  band  of  clearing  south  of  Point  Con¬ 
ception  as  a  result  of  the  subsidence  and  mixing  associated  with  the  Catalina  Fddy.  The 
stratocumulus  cloud  field  is  aligned  in  streets  revealing  the  backing  nature  of  the  surface 
wind  as  it  progresses  around  the  west  and  south  sides  of  the  eddy.  A  tongue  of  clearing 
exists  to  the  northwest,  extending  south  and  nearly  reaching  the  region  of  thin 
stratocumulus  seen  in  the  southern  extreme  of  the  image.  Another  clear  feature  is  a 
frontal  zone  indicating  the  morning  extent  of  the  diurnal  land  breeze.  As  the  drier  air 
from  the  continent  meets  the  general  flow  over  the  ocean,  upwelling  in  the  MABL  cre¬ 
ates  a  line  of  brighter  clouds  with  an  adjacent  strip  of  clouds  of  lower  optical  depth  to 
the  immediate  west.  Convergence  of  southwesterly  flow  with  the  frontal  zone  increases 
the  brightness  of  the  cloud  field  in  the  southeast  corner  of  the  image. 

The  channel  3  image  reveals  a  myriad  of  interesting  features.  Four  clear 
ship  tracks  appear  in  the  image,  with  two  more  possible  tracks  just  to  the  west  of  Point 
Conception.  A  double  track  is  found  passing  through  the  center  and  to  the  southwest. 
Note  that  these  tracks  can  still  be  traced  through  the  more  broken  stratocumulus  to  the 
southwest.  A  ship  appears  to  have  altered  course  in  one  case;  the  right-angled  track  is 
found  in  the  southeast  portion  of  the  subscene  coincident  with  cross-section  L.  Areas 
that  are  expected  to  have  smaller  cloud  droplet  size  distribution  due  to  evaporation  and 
dissipation  appear  bright,  such  as  in  the  vicinity  of  clear  areas  within  the  eddy  spiral 
band,  around  the  edges  of  dry  air  intrusion  in  the  northwest,  and  within  the  stratus  deck 
over  the  Los  Angeles  basin.  The  complexity  of  the  flow  regime  east  of  a  line  connecting 
San  Miguel,  San  Nicolas,  and  San  Clemente  Islands  due  to  the  land  breeze  and  the 
Catalina  Eddy  is  very  evident  in  the  channel  3  image.  The  land  breeze  frontal  zone  ap¬ 
pears  as  expected  from  an  examination  of  the  visual  imagery.  The  upwclled  areas  may¬ 
be  darker  as  a  result  of  larger  cloud  droplets,  while  the  brighter  region  is  consistent  with 
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Fig-  59.  San  Nicolas  Island  Rawinsonde  Sounding:  6  July  1987  1646  UTC. 


smaller  droplet  spectra  and  the  influence  of  continentally  influenced  aerosols  in  this  re¬ 
gion. 

c.  Reflectance  Cross-Sections 

Fig.  62  displays  the  reflectance  across  the  land  breeze  frontal  feature,  iden¬ 
tified  as  I  in  the  images.  Channel  1  reflectance  drops  abruptly  at  the  frontal  interface, 
from  approximately  0.85  to  0.42,  and  then  rises  past  its  original  brightness  over  the  next 
30  pixels.  Reflectance  in  channel  3  rises  from  an  average  of  0.27  to  0.35  across  the 
interface. 

Reflectance  values  for  cross-section  J  are  shown  in  Fig.  63.  This  profile 
intersects  two  small  clear  regions  within  the  Catalina  Eddy  spiral  feature.  The  rise  and 
fall  in  channel  3  reflectance  at  cloud  edges  from  decreasing  droplet  size  and  decaying 
cloud  optical  depth,  respectively,  is  obvious  in  the  channel  3  profile  and  is  similar  to  that 
observed  in  previous  case  studies.  The  rise  in  reflectance  in  both  channels  due  to  a  ship 
track  is  evident  at  the  fifth  pixel. 


Fig.  60. 


NOAA-IO  AVHRR  1524  UTC  6  July  1987  Channel  1:  Line  segments 
indicate  reflectance  cross-sections  I  through  L.  Center  coordinates 
32°35.2'N  122033.4'W. 
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NOAA-IO  AVHRR  1524  UTC  6  July  1987  Channel  3  Reflectance:  Line 
segments  indicate  reflectance  cross-sections  I  through  L.  Center  coordi¬ 
nates  32°35.2'N  122°33.4'W. 
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Fig.  62.  Reflectance  Cross-Section  I:  Sampled  west  to  east. 

Fig.  64  shows  two  relatively  short  (25  pixel)  cross-sections,  labeled  K  and 
L.  Cross-section  K  provides  reflectance  data  for  the  region  near  the  clear  intrusion  in 
the  northwest  portion  of  the  subscene,  and  intersects  a  strong  ship  track  (0.08  rise  in 
channel  3  reflectance).  The  channel  3  profile  shows  a  peak  near  the  cloud  edge  at  the 
tenth  pixel  and  a  decrease  to  the  cloud  edge  between  the  fourth  and  ninth  pixels,  con¬ 
sistent  with  other  case  studies.  Cross-section  L  bisects  another  ship  track,  and  displays 
a  0.05  rise  in  channel  3  reflectance  at  the  ship  track  position.  Both  ship  tracks  in  Fig. 
64  have  less  of  a  signature  in  channel  1  than  the  ship  tracks  seen  in  the  previous  case 
studies. 

4.  CS4  (18  July) 

a.  Synoptic  Summary 

The  west  coast  at  500  mb  is  being  influenced  by  a  deepening  wave  centered 
at  43  °N  120°W,  shown  in  Fig.  65.  The  1037  mb  subtropical  high  is  located  near  the 
same  latitude  at  147°W,  while  thermal  troughing  is  nearly  nonexistent.  The  decay  in 
subsidence  strength  in  association  with  the  trough  aloft  is  reflected  in  the  rawinsonde 
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Fig.  63.  Reflectance  Cross-Section  J:  Sampled  west  to  east;  ST  denotes  position 
of  ship  track. 


sounding  from  Vandenburg  Air  Force  Base  shown  in  Fig.  66.  The  products  of  these 
synoptic  ingredients  are  moderate  northwesterly  flow  and  scattered/ broken  cloudiness 
within  the  operating  area.  These  conditions  are  optimum  for  examining  the  reflectance 
characteristic  of  individual  cloud  elements. 
b.  A  VHRR  Subscene 

The  satellite  subscenc  images  for  CS4  are  displayed  in  Figs.  67  and  68,  and 
originate  from  the  NOAA-IO  decending  pass  crossing  35°N  at  16:04:20  UTC.  The  sub¬ 
scene  center  coordinates  are  35°16.8'N  124°12.4'W.  The  channel  I  image  reveals  a  solid 
field  of  large  cloud  cells  in  the  western  portion  of  the  subscene.  A  band  of  clouds  ap¬ 
pears  near  subscene  center,  with  small,  individual  cells  south  of  this  feature.  Closer  to 
land,  the  cloudiness  is  thinner  and  less  well-defined,  with  evidence  of  gravity  waves 
propagating  within  the  field.  The  individual  cells,  as  expected,  appear  bright  in  channel 
1  near  the  center  and  decrease  in  brightness  towards  the  cloud  edge.  The  cloud  bound¬ 
aries  are  apparent  in  channel  3  where  they  are  defined  by  high  reflectance. 

103 


>„V L> 


Fig.  64.  Reflectance  Cross-Sections  K  and  L:  Cross-section  K  sampled  west  to 
east;  L  sampled  north  to  south.  ST  denotes  position  of  ship  track. 


c.  Reflectance  Cross-Sections 

Figs.  69  and  70  display  two-dimensional  reflectance  cross-sections  within 
the  18  July  subscene,  identified  as  V  and  W,  respectively.  These  cross-sections  encom¬ 
pass  individual  cloud  cells  in  the  region  of  scattered  cloud  elements  south  of  the  subscene 
center.  Two  discrete  cells  are  within  cross-section  V,  while  cross-section  W  highlights  a 
single  element.  Channel  1  reflectance  peaks  near  the  center  of  each  element,  and  de¬ 
creases  towards  the  edge  as  EQLWC  and  optical  depth  decrease.  McKee  and  Klehr 
(1978)  shows  that  cloud  microphysical  properties  have  a  relatively  small  effect  on 
reflectance  in  visible  wavelengths  compared  to  variations  resulting  from  cloud  geometry. 
Based  on  this  study  and  the  conclusions  drawn  from  model  cloud  reflectances,  it  is  clear 
that  changes  in  cloud  thickness  and  EQLWC  are  responsible  for  most  of  the  drop  in 
reflectance  towards  the  cloud  edge  in  both  channel  1  cross-sections. 

In  channel  3,  reflectance  increases  towards  the  cloud  edge  away  from  the 
element  center  at  all  azimuths  until  at  the  edge  where  there  is  a  short  decrease.  These 
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Fig.  66.  Vandenburg  AFB  Rawinsonde  Sounding:  18  July  1987  1200  UTC. 

effects  have  been  shown  in  previous  cross-sections  and  may  result  from  the  decrease  in 
cloud  droplet  size  with  decreasing  distance  from  the  cloud  edge,  except  at  the  immediate 
edge  where  the  loss  of  optical  depth  dominates  the  microphysical  effect.  Coakley  and 
Davies  { 1986)  illustrates  that  increases  in  channel  3  reflectance  at  cloud  edges  may  be 
due  to  two  factors:  directionally  dependent  reflectance  from  cloud  sides  and  cloud-sun 
geometry,  and  variations  in  cloud  droplet  size  distribution.  The  solar  zenith  angle  at  the 
two  cross-sections  is  approximately  67°;  if  cloud  geometry  were  dominating  the  channel 
3  reflectance  variations  there  should  be  a  strong  azimuthal  dependency  of  reflectance 
circumscribing  a  cloud  element  near  the  solar  azimuth  angle  axis  (approximately  106° 
from  north).  Although  some  evidence  of  this  effect  is  visible  in  the  location  of  the 
channel  1  reflectance  maxima  in  cross-section  V,  the  effect  is  not  apparent  in  either 
cross-section  in  channel  3.  Based  on  this,  the  results  of  CS4  indicate  that  the  cloud 
droplet  size  distribution  gradient  near  the  cloud  edge  may  be  the  source  of  a  majority 
of  the  reflectance  increase  near  the  cloud  edge  in  scattered  stratocumulus/cumulus 
clouds. 
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Fig.  67.  NOAA-IO  AVHRR  1604  UTC  18  July  1987  Channel  1:  Boxes  indicate 
two-dimensional  reflectance  cross-sections  V  and  W.  Center  coordinates 
35°16.8'N  124°  1 2.4' W. 
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Fig.  68.  NOAA-IO  AVHRR  1604  UTC  18  July  1987  Channel  3 
Reflectance:  Boxes  indicate  two-dimensional  cross-sections  V  and  W. 
Center  coordinates  35°16.8'N  124°12.4'W. 


I 


108 


Fig.  69.  Tno-Dimensionai  Reflectance  Cross-Section  V:  Channel  1  (top)  and 
channel  3  reflectance  (bottom).  2-D  pixel  origin  is  to  the  northwest. 


Three  individual  cloud  elements  of  10-15  pixel  radius  are  selected  in  the  vi¬ 
cinity  of  the  aircraft  mission  track,  below  the  southern  boundary  of  the  18  July  subscene. 
This  element  size  is  selected  in  order  to  approximate  the  cloud  cell  from  which  aircraft 
measurements  are  sampled  (as  illustrated  in  the  next  section).  Pixel  values  in  these  three 
elements  are  identified  from  cloud  edge  to  cloud  center,  and  the  resulting  reflectances 
arc  averaged  over  eight  azimuth  angles  and  all  elements.  These  values  are  combined  into 
an  average  cloud  reflectance  along  a  line  segment  from  cloud  edge  to  cloud  center.  An 
examination  of  this  element  "average"  reveals  channel  1  reflectance  values  of  approxi¬ 
mately  0.25,  0.61, 0.75,  and  0.9  for  four  equal  segments  along  the  radius  line,  from  cloud 
edge  to  cloud  center.  Channel  3  reflectance  averages  0.18  at  the  cloud  edge,  increasing 
to  0.29  immediately  inside  the  cloud  edge,  with  a  near-linear  decrease  to  0.23  in  the  vi¬ 
cinity  of  the  cloud  element  center. 
d.  Aircraft  Data 
( Ij  IS  July  1987 


TABLE  12 

SELECTED  NCAR  ELECTRA  DATA  TOR  18  JULY  1987 


Time 

(LTC) 

(±20sec) 

APN-159 

Altitude 

(m) 

PMS-King 

EQLWC 

(S/'H3) 

ASASP 

cone. 

{SI  Ctrl') 

ASASP 

mean 

radius 

(^m) 

FSSP 

cone. 

(A/ cm3) 

FSSP 

mean 

radius 

(Atm) 

1453 

935 

0.22 

800 

0.085 

71 

4.0 

1454 

930 

0.29 

900 

0.090 

63 

6.0 

1454 

935 

0.16 

750 

0.090 

76 

5.3 

1517 

980 

0.16 

420 

0.095 

53 

3.1 

1517 

1050 

0.24 

405 

0.095 

75 

5.3 

1519 

1090 

0.27 

495 

0.092 

88 

6.2 

1521 

1140 

0.29 

480 

0.096 

111 

7.2 

The  flight  track  for  the  18  July  NCAR  Electra  mission  is  found  in  F  ig. 
71.  In  the  vicinity  of  the  aircraft  track,  cloud  bases  and  average  cloud  tops  were  at  al¬ 
titudes  of  900  and  1150  meters,  respectively.  Due  to  the  difference  in  time  between  air¬ 
craft  in-cloud  sampling  and  the  NOAA-IO  overpass,  a  one-to-one  match  with  a 
particular  element  is  not  attempted.  Since  the  intent  of  CS4  is  to  examine  the  reflectance 
properties  of  individual  cloud  elements,  the  aircraft  data  are  examined  in  order  to  locate 
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Fig.  71.  NCAR  Electra  Flight  Track,  18  July  1987:  (From  Kloesel  ei  ai,  1988) 


sampling  tracks  that  intersect  these  elements.  The  cloud  microphysical  data  found  in 
Table  12  is  acquired  from  the  aircraft  track  between  points  G  and  H,  uuring  an  in-cloud 
turbulence  leg  from  1452-1500  UTC  and  from  a  cloud  sounding  leg  from  1511-1523 
UTC.  In  the  first  leg,  a  small  cloud  element  is  penetrated  at  an  altitude  of  940  meters, 
with  the  aircraft  entering  the  element  at  14:52:55  UTC  and  exiting  the  other  side  at 
14:54:20  UTC.  In  the  second  leg,  the  sounding  profile  intersects  a  larger  cloud  element 
and  provides  cloud  samples  from  a  point  near  the  cloud  edge  (980  meter  altitude  at  time 
15:16:50  UTC)  to  a  region  near  the  element  center  (1140  meter  altitude  at  time  15:20:20 
UTC)  at  vertical  positions  immediately  below  the  cloud  top. 

Both  samples  indicate  decreases  in  mean  droplet  size  near  the  cloud 
edge.  Since  the  size  of  the  element  sampled  in  the  sounding  leg  is  representative  of  the 
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elements  from  which  reflectance  data  were  derived  in  the  previous  section,  the  droplet 
size  data  from  this  leg  will  be  used  for  comparison  with  averaged  channel  3  reflectances 
from  CS4.  The  short  duration  of  in-cloud  time  for  the  first  sampling  leg  and  the  rapid 
change  of  cloud  microphysical  properties  with  time  prevent  the  one-minute  averages  of 
droplet/acrosol  concentration  in  the  histograms  from  providing  useful  information  about 
the  cloud  structure.  The  first  leg  reveals  the  expected  decrease  in  cloud  droplet  size, 
concentration,  and  EQLWC  near  the  edge  of  the  cloud,  while  the  aerosol  parameters 
remain  relatively  unchanged.  Histograms  for  the  cloud  edge  and  the  cloud  interior  in  the 
second  sampling  leg  are  shown  in  Fig.  72.  The  changes  in  cloud  droplet  size  distribution 
across  the  cloud  element  are  clearly  evident.  As  expected,  the  distribution  of  aerosols 
at  both  locations  are  nearly  identical,  while  the  droplet  distrbution  shows  an  unmistak¬ 
able  shift  toward  smaller  droplet  sizes  near  the  cloud  edge.  The  absence  of  large  cloud 
droplets  and  the  concentration  of  aerosol  particles  indicate  that  the  sampled  cloud  ele¬ 
ments  are  within  continentally  influenced  air. 

5.  Reflectance  Analysis 

In  order  to  correlate  the  cloud  droplet  size  data  from  the  NCAR  Electra  with 
reflectance  data  from  applicable  NOAA-9/10  subscenes,  the  aircraft  samples  listed  in 
Tables  9  through  1 1  are  examined  to  determine  the  samples  that  coincide  with  the  ap¬ 
proximate  location  of  reflectance  cross-section  data.  When  a  match  is  identified,  a  ten- 
pixel  average  of  reflectance  centered  near  the  coincident  pixel  is  performed.  This  average 
reflectance  value  is  utilized  to  make  comparisons  with  the  cloud  microphysical  data.  The 
limits  of  uncertainty  represent  a  95%  confidence  (within  two  standard  deviations)  in  the 
value  of  reflectance  along  the  pixel  average,  and  in  the  droplet  size  determination.  The 
uncertainty  of  droplet  size  values  are  a  function  of  the  variability  and  noise  contained 
in  the  aircraft  data.  In  the  case  of  the  18  July  case  study,  the  cloud  element  averages 
identified  in  the  previous  section  are  correlated  to  the  aircraft  sounding  data  contained 
in  Table  12. 

The  results  of  the  intercomparisons  are  displayed  in  Fig.  73.  As  expected,  there 
is  little  correlation  between  reflectance  and  cloud  droplet  size  in  channel  1  since  the 
magnitude  of  reflectance  in  this  portion  of  the  electromagnetic  spectrum  is  heavily  de¬ 
pendent  upon  cloud  thickness  and  EQLWC,  as  well  as  the  droplet  size  distribution.  The 
four  data  points  between  3  and  7jxm  represent  individual  cloud  element  averages  from 
CS4.  The  data  show  an  increase  in  reflectance  as  larger  droplets  are  encountered  in  the 
interior  of  the  cloud  element  where  greater  EQLWC  and  cloud  thickness  exists.  The 


Fig.  72. 
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Size  Distribution  Histograms  for  18  July  1987:  ASASP  aerosol  size  dis¬ 
tribution  (top)  and  FSSP  cloud  droplet  size  distribution  (bottom). 
Dashed  lines  indicate  samples  near  cloud  element  interior,  solid  lines  rep¬ 
resent  the  cloud  edge  regime. 
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_ FSSP  MEAN  DROP  RADIUS  (MICRON! _ 

Reflectance  verses  FSSP  Mean  Droplet  Radius:  Channel  l  (top)  and 
channel  3  (bottom).  Brackets  indicate  95  percent  confidence  intervals  for 
reflectance  and  mean  radius. 
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channel  3  reflectance  displays  a  near-linear  correlation  with  cloud  droplet  size.  The 
outlying  data  point  at  3/xm  represents  the  CS4  cloud  element  edge,  where  the  decrease 
in  cloud  thickness  and  its  associated  effect  on  the  optical  depth  is  dominating  the 
microphysical  effect  on  reflectance.  A  review  of  the  aircraft  data  indicates  that  the  cloud 
thickness  for  this  element  edge  sample  is  less  than  SO  meters,  and  the  LQLWC  is  almost 
half  of  that  found  in  the  cloud  interior  (0.16  g/m*  verses  0.29  g/m1  ). 

Fig.  74  relates  channel  1  reflectance  to  cloud  thickness.  Cloud  thickness  values 
are  estimated  from  APN-159  radar  altimeter  and  PMS-King  LQLWC  data.  The  esti¬ 
mations  may  vary  from  actual  cloud  thicknesses  by  several  tens  of  meters  since  aircraft 
measurements  must  be  interpolated  horizontally  in  certain  aircraft  mission  data.  The 
data  show  increasing  channel  1  reflectance  with  increasing  cloud  thickness,  which  is 
consistent  with  the  dependence  on  optical  depth  at  this  wavelength.  Although  channel 
1  reflectance  should  vary  as  a  function  of  droplet  radius  with  constant  cloud  thickness 
and  LQLWC,  there  are  insufficient  data  to  observe  this  variation. 

Fig.  75  compares  the  channel  3  results  with  the  model  cloud  reflectance  data 
discussed  and  illustrated  in  Chapter  II.  Noonkester  (1984)  showed  that  the  cloud 
droplet  distribution  modal  radius  is,  on  the  average,  approximately  lpim  larger  than  the 
mean  droplet  radius;  this  is  consistent  with  the  NCAR  Electra  data  utilized  in  this  thesis. 
The  difference  between  droplet  distribution  modal  radius  (used  in  model  reflectances) 
and  droplet  distribution  mean  radius  (obtained  from  the  aircraft  data)  adds  a  small  de¬ 
gree  of  uncertainty  to  the  data  in  Fig.  75,  and  is  considered  in  the  definition  of  the  region 
of  uncertainty.  The  model  case  with  a  cloud  thickness  of  250  meters  and  EQLWC  of 
0.4  g!ml  is  the  most  representative  of  the  conditions  specified  in  the  model  reflectance 
data  in  Chapter  II.  For  this  comparison,  three  additional  model  reflectances  were  gen¬ 
erated  for  a  droplet  modal  radius  of  12um  in  order  to  compare  delta-Eddington  ap¬ 
proximations  with  the  aircraft  and  reflectance  data  at  larger  droplet  sizes.  As  seen  in 
Fig.  75,  a  majority  of  the  region  of  uncertainty  falls  within  the  region  defined  by  the 
three  model  cloud  droplet  distributions.  The  shift  toward  model  distribution  D1  at  large 
droplet  sizes  and  the  shift  toward  model  distribution  D2  at  the  smaller  droplet  sizes  is 
consistent  with  the  shift  in  distribution  shape  measured  by  the  aircraft  and  illustrated  in 
Figs.  33  and  57.  Based  upon  the  trend  observed  at  large  droplet  radii,  the  relationship 
between  droplet  size  and  channel  3  reflectance  will  continue  to  asymptotically  approach 
a  reflectance  value  of  zero  for  droplet  radii  larger  than  the  largest  mean  droplet  radius 
measured  in  the  case  studies. 


CLOUD  THICKNESS  (M) 


Fig.  74.  Channel  1  Reflectance  verses  Cloud  Thickness:  Cloud  thicknesses  esti¬ 
mated  from  APN-159  radar  altimeter  and  PMS-King  EQLWC  data. 
Brackets  indicate  95  percent  confidence  intervals  for  reflectance  and  cloud 
thickness. 
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RADIUS  (MICROMETERS) 


Channel  3  Reflectance  verses  Droplet  Radius:  Dashed  lines  indicate 
model  cloud  reflectance  for  droplet  distributions  Dl,  D2,  and  D3;  solid 
lines  indicate  95  percent  confidence  interval  for  data  points. 


riy«HWunuournjfW7vnv Jnoronewv  w  v. v  vjv.va.va  wnuu  v*  v*  uv  ™  «  ««  m.  ir* 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 
A.  CONCLUSIONS 

The  relationship  between  solar  reflectance  and  cloud  droplet  size  distribution  has 
been  illustrated  by  comparing  AVIIRR  satellite  data  and  aircraft  measurements.  The 
expected  relationship  of  higher  reflectances  from  smaller  cloud  droplet  size  spectra  is 
confirmed  by  the  four  case  studies  presented  in  this  thesis.  In  the  case  of 
continental  marine  air  mass  differences,  higher  channel  3  reflectivities  are  associated 
with  higher  concentrations  of  aerosols,  higher  concentrations  of  cloud  droplets,  and  a 
shift  towards  a  smaller  mean  cloud  droplet  radius.  In  the  case  of  individual  cloud  ele¬ 
ments,  higher  channel  3  reflectivities  are  associated  with  a  strong  shift  towards  a  smaller 
mean  droplet  radius,  even  though  EQLWC  and  droplet  concentration  decrease  and  the 
aerosol  concentration  remains  relatively  unchanged.  The  strength  of  channel  1 
reflectance  provides  a  measure  of  optical  depth,  ensuring  that  it  is  sufficiently  "thick"  to 
allow  deductions  of  cloud  properties  from  channel  3  measurements.  The  strong  de¬ 
pendence  of  channel  3  reflectance  on  cloud  droplet  size  distribution  allows  inference  of 
cloud  composition  characteristics  from  satellite  data.  This  dependence  is  closely  ap¬ 
proximated  by  model  cloud  reflectances  based  on  Mie  theory  and  the  delta-Eddington 
approximation.  Assuming  that  the  MABL  is  well-mixed,  the  satellite  reflectance  data 
provide  an  indication  of  the  concentration  and  nature  of  the  aerosol  distribution  below 
the  cloud  layer. 

Regions  of  continental  air  mass  influence  over  the  ocean  may  be  generated  by 
synoptically  forced  offshore  flow  over  the  western  coast  of  the  U.S.  and  Canada,  and 
may  be  advected  with  the  wind  field.  These  regions  may  vary  in  their  concentration  of 
continental  aerosols  and  associated  CCN  (and  mean  cloud  droplet  size)  from  the 
strength  of  the  offshore  flow,  as  well  as  the  duration  of  time  that  this  flow  persists  prior 
to  it  being  replaced  by  the  climatological  northwesterly  flow  from  marine  sources.  The 
continental,  marine  air  mass  interface  identified  in  the  satellite  data  seems  to  be  advected 
at  a  rate  slightly  slower  than  the  wind  forcing,  suggesting  that  some  mixing  of  aerosols 
persists  across  the  interface.  The  1550  per  cm 3  aerosol  concentration  measured  at  a 
point  approximately  300  n  mi  from  the  coastline  in  CS2  indicates  that  significant  con¬ 
centrations  of  continental  aerosol  have  the  potential  to  be  advected  far  offshore  and  af¬ 
fect  the  reflectance  of  stratocumulus  cloud  fields  over  thousands  of  square  miles. 
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The  stratocumulus  cloud  fields  within  the  continental  influence  may  appear  more 
continuous  and  less  broken  than  clouds  within  the  marine  influence  under  similar  con¬ 
ditions.  In  several  cases,  continental,  marine  interfaces  identified  in  channel  3  reflectance 
imagery  can  be  located  utilizing  the  variations  observed  in  channel  1. 

Although  model  cloud  reflectance  calculations  suggest  that  continental!)'  influenced 
clouds  should  increase  in  brightness  in  channel  1  as  well  as  channel  3,  these  regions  ap¬ 
pear  less  bright  than  the  marine  stratocumulus  in  the  visible.  Measurements  of  cloud 
top  EQLWC  in  both  air  mass  regimes  indicate  relatively  constant  values  (near  0.31 
gjm 3)  which  are  consistent  with  those  observed  by  Noonkcstcr  (198-1).  Only  cloud-edge 
measurements  give  significantly  lower  levels  of  EQLWC.  It  appears  that  changes  in 
stratocumulus  cloud  thickness  and  horizontal  cloud  coverage  occur  in  cloud  fields  mod¬ 
ified  by  the  addition  of  continental  CCN,  so  that  the  scattering  optical  depth  is  reduced 
in  the  continental  regime.  This  reduces  the  reflectance  in  channel  1  but  does  not  afTect 
channel  3  reflectance  for  stratocumulus  clouds  of  sufficient  thickness  (>  100  meters). 
The  model  reflectance  calculations  seem  to  hold  true  for  ship  effluent  tracks,  where  in¬ 
creased  reflectivity  in  both  satellite  channels  is  observed. 

Cloud  elements  surrounded  by  sub-saturated  air  exhibit  a  characteristic  signature  in 
channel  3  reflectance.  This  signature  is  consistent  with  horizontally  vary  ing  droplet  size 
distribution  and  optical  depth.  The  channel  3  reflectance  increases  from  the  cloud  cle¬ 
ment  center  to  the  cloud  edge  where  the  greatest  entrainment  of  ambient  air  is  taking 
place.  The  reflectance  decreases  at  the  immediate  cloud  edge  where  the  loss  of  optical 
depth  dominates  the  inicrophysical  effect.  This  signature  appears  to  be  of  a  greater  in¬ 
fluence  on  channel  3  reflectance  than  variations  in  cloud  and  sun-cioud-satellite  geom¬ 
etry',  and  is  even  observed  in  individual  convective  turrets  within  a  stratocumulus  cloud 
field.  The  effect  is  also  observed  at  cloud  edges  other  than  those  resulting  from  the 
convective  circulation  within  a  cloud  element,  such  as  the  clear  areas  due  to  tuibulcnt 
mixing  associated  with  the  Catalina  Eddy  in  CS3.  Channel  l  reflectance  measurements 
in  cloud  elements  show  a  decrease  in  reflectance  towards  the  cloud  edge  as  a  result  of 
decreasing  optical  depth. 

Model  cloud  reflectance  calculations  suggest  that  higher  solar  reflectances  occur 
with  droplet  size  distributions  that  have  a  strong  peak  at  the  modal  radius,  while  lower 
reflectance  results  from  a  broader  and  less  sttongly  peaked  droplet  size  distribution.  The 
trend  observed  in  the  reflectance/aircraft  data  support  this  dependence.  The  mean 
droplet  size  measurements  in  marine  environments  correspond  to  the  larger  droplet  sizes 


which  are  shifted  in  reflectance  towards  the  values  represented  by  the  broad  model 
droplet  size  distribution;  while  continentally  inlluenced  values  arc  shifted  toward  the 
strongly  peaked  model  droplet  size  distribution  reflectances.  The  droplet  spectra 
histograms  reveal  the  change  between  environments  and  the  consistency  with  the  mod¬ 
eled  expectations. 

The  relationship  between  channel  3  reflectance  and  mean  cloud  droplet  radius  illus¬ 
trated  in  this  thesis  indicates  the  potential  to  make  accurate  estimations  of  the  micro¬ 
physical  properties  of  the  MABL  based  upon  satellite  data.  Even  with  the  limited 
sample  contained  herein,  mean  cloud  droplet  sizes  can  be  estimated  within  an  approxi¬ 
mate  range  of  3nm  with  95%  confidence  for  mean  cloud  droplet  sizes  less  than  lOyum  in 
stratocumulus  clouds  of  sufficient  optical  depth.  At  larger  mean  radii,  the  small  change 
in  channel  3  reflectance  with  varying  droplet  size  make  the  estimation  more  difficult. 
The  precision  of  this  estimation  should  improve  as  more  reflectance/droplet  size  corre¬ 
lations  are  made  in  future  research. 

Observations  of  reflectance  from  channels  1  and  3  have  the  potential  to  reveal  the 
location,  extent,  and  the  horizontal  microphysical  variability  of  sub-synoptic  scale  phe¬ 
nomena  such  as  land-sea  breeze  circulations,  topographically  induced  features  such  as 
the  Catalina  Eddy,  gravity  waves,  and  sources  of  natural  and  industrial  pollution.  Such 
a  capability  may  provide  a  tool  to  ascertain  the  impact  of  events  such  as  forest  fires, 
industrial  accidents,  or  volcanic  activity  on  the  environment. 

B.  RECOMMENDATIONS 

Future  observations  of  marine  stratocumulus  clouds,  such  as  those  which  will  occur 
in  conjunction  with  the  next  FIRE  marine  stratocumulus  IFO  tentatively  scheduled  for 
July  1990  at  an  undetermined  site,  should  include  plans  to  vector  aircraft  across  strong 
continental  marine  air  mass  boundaries  on  short  notice  by  experiment  team  members 
collecting  real-time  AVHRR  data.  Ideally,  the  aircraft  would  perform  a  cloud  sounding 
leg  of  100+  nautical  miles  centered  on  the  air  mass  interface,  with  the  sounding  leg  be¬ 
ing  coincident  with  the  satellite  overpass  subsequent  to  the  overpass  used  to  localize  the 
region  of  interest.  Instructions  concerning  the  location  of  applicable  sampling  locations 
could  be  relayed  to  aircraft  enroute  to  its  geographic  objective.  Likewise,  samples 
through  individual  cloud  elements  should  be  performed  as  close  as  possible  to  satellite 
overpass  times.  This  would  allow  precise  correlation  between  measured  reflectance  and 
cloud  elements,  and  would  provide  the  opportunity  for  measurements  to  be  performed 
in  local,  short-duration  phenomena  such  as  ship  effluent  tracks.  This  would  also  avoid 


the  "hit-or-miss"  strategy  of  sampling  specific  features  within  the  stratocumulus  cloud 
field. 

The  connection  between  channel  3  reflectance  and  MABL  turbulence  needs  to  be 
examined.  Entity  type  entrainment  mixing  (ETEM),  discussed  by  Telford  and  Wagner 
(1981)  and  Telford  et  al.  (1984),  may  describe  the  interaction  between  the  stratocumulus 
cloud  top  and  the  subsidence  inversion.  ETEM  in  a  stratocumulus  cloud  regime  changes 
little  in  the  horizontal  in  the  absence  of  topographic  influences.  The  entrainment  results 
in  a  gradient  of  cloud  droplet  sizes  at  the  cloud  top,  with  the  depth  of  this  gradient  de¬ 
pendent  upon  the  strength  of  turbulent  energy  within  the  MABL.  This  droplet  size 
gradient  is  within  the  depth  of  the  cloud  acting  as  the  source  of  channel  3  reflectance. 
It  is  possible  that  changes  in  cloud  top  entrainment  due  to  turbulence  variations  caused 
by  changes  in  sea-surface  temperature,  MABL  wind  speed,  and  synoptic  forcing  could 
alfect  channel  3  reflectance  sufficiently  enough  to  reduce  the  precision  of  cloud  droplet 
size  estimation  techniques.  It  has  been  shown  that  entrainment  near  the  edge  of  indi¬ 
vidual  cloud  elements  has  a  significant  effect  on  channel  3  reflectance.  Accounting  for 
the  effect  of  turbulence  on  solar  reflectance  will  allow  estimations  of  cloud  droplet  size 
distribution  to  be  more  closely  correlated  to  the  aerosol  environment  within  the  MABL 
since  the  CCN  composition  would  dominate  the  factors  affecting  cloud  droplet  size  dis¬ 
tribution. 

Variations  in  reflectance  due  to  anisotropy  need  to  be  described  in  greater  detail. 
As  mentioned  in  Chapter  II,  the  anisotropic  reflectance  factors  provided  by  Taylor  and 
Stowe  (1984)  represent  spectrally  integrated  values  for  wavelengths  between  OZ/nn  and 
50£tm,  are  based  on  a  finite  number  of  NIMBUS-7  data,  and  represent  an  average  value 
for  a  discrete  angular  bin.  It  is  likely  that  the  anisotropic  reflectance  factors  vary  as  a 
function  of  wavelength.  Regardless  of  the  quality  of  any  smoothing  routine  designed  to 
approximate  the  effect  of  anisotropy  over  the  entire  domain  of  reflective  surfaces,  satel¬ 
lite  and  solar  zenith  angles,  and  relative  azimuth  angles,  it  will  only  be  as  good  as  the 
accuracy  of  the  original  bin  values.  As  the  availability  of  anisotropic  reflectance  data 
grows,  the  revised  measurements  should  be  incorporated  into  any  scheme  dependent 
upon  the  conversion  of  directional  reflectance  to  isotropic  reflectance. 

The  effect  of  changing  cloud  reflectance  due  to  spatial  variations  in  aerosol  con¬ 
centration  and  composition  on  the  global  climate  must  be  incorporated  into  climate 
models.  As  discussed  by  Twomey  et  al.  (1984),  slight  variations  in  cloud-related  pa¬ 
rameters  in  global  climate  models  results  in  a  considerable  climatic  response.  Failure  to 


consider  clouds  as  reflectance  "variables"  may  lead  to  an  underestimate  of  global  albedo. 
The  simplicity  of  the  technique  utilized  in  this  thesis  makes  it  ideal  for  performinu 
climatological  studies  through  the  use  of  archived  AVHRR  data.  Through  a  study  such 
as  this,  temporal  trends  in  cloud  reflectance  due  to  microphysical  variations  could  be 
identified.  This  potential  may  be  especially  valuable  in  describing  the  influence  of  D.V1S 
on  cloud  albedo  over  remote  oceanic  regions  where  in-situ  measurements  are  impracti¬ 
cal. 

Model  reflectance  data  should  be  assimilated  in  preparation  for  the  new  1.6um  im¬ 
aging  channel  to  be  included  with  the  AVHRR  instrument  on  future  satellites.  As  sug¬ 
gested  by  Wetzel  and  Vondcr  Ilaar  (19S6).  the  implementation  of  this  new  sensor 
coupled  with  the  potential  for  testing  and  utilizing  more  sophisticated  sensors  on  De¬ 
fense  Meteorological  Satellite  Program  (DMSP)  and  GOES  platforms  provides  impetus 
for  increasing  locus  on  the  applications  of  N1R  satellite  imagery  in  determining  cloud 
micropin sical  structure.  The  ability  to  cross-correlate  NIR  reflectance  from  two  sepa¬ 
rate  wavelength  bands  will  enable  more  piecise  determinations  of  cloud  microphysical 
characteristics. 
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APPENDIX  COMPOSITE  IMAGES 

The  composite  RGB  images  [or  each  of  the  nine  case  study  subscenes  are  shown  in 
Fig.  76  through  84.  The  line  segments  indicate  reflectance  cross-sections  corresponding 
to  those  which  appear  in  the  monochromatic  imagery  in  Chapter  IV.  Chapter  111  ex¬ 
plains  the  interpretation  of  these  images. 
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Fig.  80.  NOAA-9  AVHRR  2321  1)7 C  13  July  1987  Composite  Image 
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82.  NOAA-9  AYHRR  2310  UTC  14  July  1987  Composite  linage 
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